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Abstract 


An  analytical  model  is  developed  to  describe  the 
response  of  the  "joggle-lap"  joint  to  both  tensile  and 
bending  loads.  The  model  consists  of  a  non-linear  beam 
analysis  which  calculates  stress  profiles  through  the 
adherent  thickness.  A  plane  stress  finite-element  model 
was  incorporated  into  the  analysis  to  correctly  determine 
the  stress  field  in  the  adhesive  zone  where  it  was  shown 
that  beam  analysis  was  less  accurate.  Elastic  response  of 
the  "joggle-lap"  joint  due  to  tensile  loads  was  verified 
through  experimental  testing  and  ultimate  loads  were 
accurately  predicted  within  experimental  error.  Maximum 
adherent  flexural  stress  was  found  to  determine  joint 
failure.  A  parametric  study  was  undertaken  by  using  the 
verified  analytical  model  and  the  results  were  recorded  as 
a  series  of  design  curves. 
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I .  Introduction 


Recent  government  regulations  for  increased 
gasoline  mileage  requirements  have  induced  automobile 
manufacturers  to  seek  light  weight  replacement  material 
systems  for  existing  metal  parts.  Since  the  automotive 
industry  is  a  high  volume  operation,  sheet  molding  compound 
(SMC)  parts  offer  a  feasible  answer  to  the  problem.  The 
SMC  molding  time  of  from  1  to  3  min/piece  depending  on  the 
size  and  thickness  of  the  part  is  compatible  with  auto¬ 
motive  assembly  line  production. 

International  Harvester  et  al  are  currently 
employing  SMC  molded  body  components  on  their  vehicles  to 
replace  former  sheet  metal  parts.  This  new  direction  has 
brought  with  it  several  problems ,  one  of  which  is  the 
design  of  adhesive  joints.  The  joint  must  accommodate 
high  rate  fabrication  techniques  and  provide  optimum 
strength  and  durability.  In  addition,  the  joint  must 
satisfy  certain  cosmetic  requirements  such  as  adjacent  , 
flush  edges.  With  these  criteria  in  mind,  the  "joggle- 
lap"  joint  has  been  chosen  for  detailed  study  and  analysis. 
This  joint  configuration  is  shown  in  Figure  1.  Since  a 
joint  of  this  type  experiences  a  variety  of  loading 
conditions  in  practice,  it  was  decided  to  model  the  joint 
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FIGURE  1:  THE  JOGGLE  LAP  JOINT  SUBJECT  TO  TENSILE  AND  BENDING  LOADS' 
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in  pure  tension  and  pure  bending.  By  superposition,  it  is 
apparent  that  any  combination  of  these  two  loading  conditions 
may  then  be  constructed. 

This  work  focuses  on  the  development  of  an  analytical 
model  to  describe  the  behavior  of  the  "joggle-lap"  joint 
due  to  both  tensile  and  bending  loading  conditions.  The 
first  section  utilizes  small  deflection  beam  theory  for 
both  straight  and  curved  beam  elements  to  obtain  a  solution 
for  the  displacement  and  stress  fields  of  the  joint. 

Included  in  this  analysis  is  the  derivation  of  the  governing 
differential  equations  for  the  deflection  of  the  curved 
beam. 

The  second  section  utilizes  a  finite-element  model 
to  reveal  localized  stress  concentrations  in  the  adhesive 
zone.  Boundary  conditions  for  the  finite  element  model  are 
obtained  from  a  transformation  of  stresses  in  the  deformed 
geometry  to  equivalent  stresses  in  the  undeformed  geometry. 
This  transformation  of  stresses  is  performed  via  a  computer 
routine  for  ease  of  calculation. 

Finally,  experimental  verification  of  the  analytical 
predictions  is  reported  along  with  a  description  of  testing 
procedures.  The  maximum  flexural  stress  is  shown  to  correlate 
strength  data  and  failure  analysis.  Also,  the  microstructure 
of  the  joint  was  examined  as  a  possible  explanation  of  the 


failure  mode. 


II-  Background 


A.  Adherent  Materials 

The  adherents  of  the  proposed  "joggle-lap"  joint 
were  composed  of  a  random-fiber  composite  known  as  SMC-25. 
SMC  is  defined  as  a  sheet  molding  compound  that  contains 
reinforcements  with  an  average  fiber  length  of  approximately 
1  inch  (2.54  cm)  with  random  orientation  in  the  plane.  The 
number  25  indicates  that  the  composite  is  25  percent  glass 
fibers  by  weight.  The  major  constituents  of  SMC  are  E-glass 
fibers  and  a  styrenated  polyester  resin  in  the  form  of  a 
paste.  It  is  quite  common  to  use  mineral  fillers  during 
the  manufacture  of  the  paste  to  facilitate  flow  when  molding 
or  to  obtain  certain  characteristics  from  the  molded  part 
such  as  a  high  resistance  to  flame  or  increased  stiffness. 
Another  prime  reason  for  using  fillers  is  the  fact  that  they 
are  much  cheaper  than  the  polyester  resin  itself  and  thus 
reduce  the  cost  of  materials.  At  times,  chemical  additives 
may  also  be  introduced  into  the  paste  to  serve  as  catalysts 
during  the  molding  cycle. 

The  process  of  SMC  manufacturing  is  a  highly 
innovative  one  which  is  completely  automated.  Figure  2 
(taken  from  Owens/Corning  Fiberglass  SMC  Review)  depicts 
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FIGURE  2:  SMC-R  MACHINE 
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a  typical  process  currently  in  use  by  a  competitive  supplier 
*  * 

of  SMC.  The  first  step  of  the  procedure  is  to  distribute 
the  resin  onto  a  polyethylene  carrier  film  as  shown.  Con¬ 
tinuous  glass  fibers  are  then  chopped  into  lengths  of  less 
than  three  inches  and  distributed  in  a  random  fashion  on 
the  wetted  film.  A  second  layer  of  resin-coated  polyethylene 
film  serves  as  a  top  layer  to  the  sandwich-like  sheet. 

Several  rows  of  rollers  act  to  insure  that  the  glass  fibers 
are  fully  impregnated  with  the  polyester  resin  thus  yielding 
consistency  in  moldability  of  the  SMC.  Finally  the  product 
is  directed  to  a  take-up  roll  for  ease  of  handling  during 
shipping  and  storage . 

SMC  is  usually  placed  in  a  constant  temperature  room 
while  storing  to  allow  maturation  to  take  place.  Maturation 
is  nothing  more  than  allowing  the  SMC  to  increase 
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in  viscosity  to  enhance  relative  ease  of  handling  of  the 
sheet.  Maturing  the  SMC  sheet  for  extended  periods  of 
time  greatly  reduces  the  flow  characteristics  of  the 
product  while  molding.  Recommended  shelf-life  for  SMC 
stored  at  10-15°  C  is  about  2  weeks,  however  in  general 
it  may  often  be  used  up  to  2  months  after  the  date  of  its 
manufacture . 

Once  the  SMC  sheet  has  reached  maturity,  it  is 
ready  for  molding.  Upon  removing  the  protective  poly¬ 
ethylene  film,  the  molding  compound  is  cut  to  size  and 
strategically  placed  in  the  mold.  This  procedure  is 
known  as  charging  the  mold.  The  so-called  strategic 
locations  of  the  mold  are  those  positions  that  allow  the 
SMC  to  flow  to  all  parts  of  the  mold  and  maintain  uni¬ 
form  part  thickness.  To  date  these  locations  have  been 
determined  by  trial  and  error  coupled  with  experience. 

Compression  molding  combines  both  temperature  and 
pressure  to  induce  an  exothermic  reaction  which  serves  to 
cure  the  part  in  the  mold.  Figure  3  (taken  from  ref  [3]) 
is  an  example  of  a  typical  curing  cycle  showing  the 
temperature  of  the  part  as  a  function  of  time .  It  should 
be  noted  that  platen  temperatures  of  200°  C  are  usually 
sufficient  for  SMC  molding  and  may  be  achieved  with  super¬ 
heated  steam.  Another  important  fact  seen  from  the  figure 
is  the  overall  cure  time .  Average  cure  times  are  generally 
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1-3  minutes  (depending  upon  the  thickness  of  the  part) 
which  lends  itself  to  production  line  applications 
inherent  in  the  automotive  industry.  Figure  4  (taken 
from  ref  13])  shows  the  effect  of  pressure  upon  a  typical 
cure  cycle.  Note  that  the  peak  pressure  and  maximum 
temperature  correspond  to  the  initiation  of  the  exothermic 
reaction.  The  key  to  successful  molding  is  to  acquire 
fine  control  of  the  application  of  pressure  to  the  cure 
cycle . 

The  main  feature  of  SMC  is  the  ability  of  the  glass 
fibers  to  flow  with  the  paste  during  the  molding  process. 
Since  the  fibers  are  transported  to  all  parts  of  the  mold, 
it  is  possible  to  produce  a  geometrically  complicated  part 
with  quasi-constant  mechanical  properties.  It  has  been 
shown  by  Pipes  and  Taggart  [ref  5] ,  that  in  areas  of 
intensified  flow,  the  fibers  tend  to  align  themselves  with 
the  direction  of  flow  and  thus  produce  areas  of  varying 
mechanical  properties.  It  is  therefore  beneficial  to  under¬ 
stand  the  flow  characteristics  within  the  mold  to  produce 
a  part  with  controlled  and/or  uniform  mechanical  properties. 
Taggart  et  al  have  determined  the  properties  of  SMC-25  to 
be  those  found  in  Table  1.  Some  scattering  in  the  data 
was  reported  due  to  the  inherent  local  variations  in  the 
material.  To  determine  the  normal  modulus  (modulus  normal 
to  the  plane  of  the  fibers) ,  the  relationship  shown  may  be 


FIGURE  4:  PRESSURE  VARIATION  OF  A  TYPICAL  CURE 
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Table  1 

Properties  of  SMC-25 


Tension 


gtension 

L 

GPa  (Msi) 

14.48 

(2.1) 

^tension 

.3 

.tension 

°ult 

MPa  (ksi) 

90 

(13.1) 

tension 

£ult 

(y  in/in) 

11,400 

Compression 

gCompression 

L 

GPa  (Msi) 

12.41 

(1.8) 

compression 

.28 

compression 

°ult 

MPa  (ksi) 

204 

(29.6) 

compression 

ult 

(y  in/in) 

20,600 

-11- 


used.  This  relationship  resembles  the  well-known  rule  of 
mixtures  for  continuous  fibrous  composites. 


where  En  =  normal  modulus  of  elasticity  of  the  composite 
vf  =  volume  fraction  of  fiber 
vm  =  volume  fraction  of  matrix 
Eg  =  modulus  of  glass  fiber 
E^,  =  modulus  of  matrix 


Table  2  provides  the  needed  data  for  determining  the 
normal  modulus  of  elasticity.  By  definition,  SMC  is  composed 


Table  2' 


Modulus  of  Elasticity 
(106  psi) 
Specific  gravity 


polyester  resin  E-glass  fiber 


1.28 


2.54 


of  25%  fiber  by  weight.  Utilizing  the  equation  written  below 


vf  +  vm  = 
t  m 


1 


(2) 


allows  one  to  solve  for  vf  where  vm  may  be  rewritten  as 


,  Sf  ,  .  ,  Wm  , 

vm  *  vf  ‘-S— ]  ‘-WT3 

m  t 


=  specific  gravity  of  fiber 
=  specific  gravity  of  matrix 
=  weight  fraction  of  matrix 
=  weight  fraction  of  fiber 


Making  the  appropriate  substitutions,  Eq.  (2)  becomes 


.75  r  2.54 
.25  L  1.28 


[  v-^-6-1  v-  +  v-  -  1 


Thus  the  corresponding  volume  fraction  of  fiber  and  matrix 
are  .14  and  .86  respectively.  From  Eq.  (1)  the  value  of  E 

g 

is  now  calculated  to  be  0.58  x  10  psi. 


LVinson  and  Chou 


B.  Adhesive  Materials 

The  adhesive  system  chosen  for  the  "joggle-lap" 
joint  was  developed  by  the  Adhesives  Division  of  Goodyear 
Chemicals.  The  Pliogrip  6000  series  is  a  general  purpose 
structural  adhesive  with  a  polyurethane  base.  Currently 
available  as  a  two-part  system,  Pliogrip  6000  exhibits 
both  high  flexibility  and  resilience.  With  the  proper 
selection  of  curatives,  the  working  time  of  the  adhesive 
may  be  accurately  controlled  between  1-6  minutes. 

In  order  to  utilize  this  adhesive  system  only 
minimal  surface  preparation  is  necessary.  The  two  surfaces 
to  be  bonded  are  prepared  with  a  plastic  wash  primer 
(Pliogrip  6033/6034  Wash  Primer)  that  is  applied  with  a 
cloth.  No  sand  blasting  or  surface  stripping  is  necessary. 
To  maintain  reliably  bonded  parts,  Pliogrip  6000  must  be 
mixed  at  a  precise  ratio  of  4  parts  resin  to  1  part  curative 
by  weight  or  volume.  Deviations  from  this  standard  will 
yield  resin-rich  areas  of  uncured  adhesive.  The  actual 
mixing  of  the  two  components  must  be  carried  out  without 
the  introduction  of  air  into  the  system,  thus  the  need  for 
specialized  equipment.  Without  this  precaution,  entrapped 
air  bubbles  in  the  cured  adhesive  would  yield  voids  and 
greatly  affect  the  performance  of  the  bond.  Curing  this 
adhesive  system  can  be  accomplished  at  room  temperature, 
however  the  use  of  heated  fixtures  will  reduce  cure  times. 
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Recommended  clamping  pressures  of  heated  fixtures  range  from 
20  to  40  psi. 

An  important  criterion  in  the  design  of  bonded  joints 
is  that  of  the  adhesive  thickness.  It  has  been  shown  that 
adhesive  properties  vary  inversely  with  adhesive  thickness. 
Thus  the  bulk  properties  of  the  adhesive  are  distinctly 
different  from  those  in  the  film  state.  So  the  question  is 
posed  as  to  the  optimum  bond  thickness  as  a  function  of  shear 
strength.  Figure  5  (taken  from  Pliogrip  technical  data, 
Goodyear  Adhesives)  shows  the  effect  of  glue  line  thickness 
on  bond  joint  strength.  A  bond  line  thickness  of  0.030 
inches  was  chosen  as  optimal  even  though  thicknesses  less 
than  0.030  inches  yield  greater  bond  strengths.  It  was  felt 
that  bonding  thicknesses  less  than  0.030  inches  are  not 
capable  of  being  fabricated  with  consistency  under  production 
operations.  (i.e.  molded  FRP  parts  will  inherently  not  fit 
together  with  reliable  precision) . 

To  achieve  uniform  bond  lines ,  one  of  two  procedures 
is  generally  used.  Adherents  may  have  a  small  raised  button 
of  0.030  inches  in  thickness  which  acts  as  a  spacer  for  the 
joint  to  insure  a  uniform  bond.  Another  procedure  is  to 
introduce  small  glass  spheres  (0.030  inches  dia.)  directly 
into  the  adhesive  to  achieve  similar  spacing.  The  effect  of 
these  spheres  on  joint  strength  has  not  been  determined  but 
it  is  argued  that  the  variation  from  the  norm  is  negligible. 


FIGURE  5:  EFFECT  OF  GLUE  LINE  THICKNESS 
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III.  Methods  of  Analysis 


A.  Tensile  Loading 
a.  Beam  Model 

Recently,  Adkins  [ref  2]  investigated  the  response 
of  a  scarf  joint  to  simple  tensile  loadings.  It  was  found 
that  the  scarf  joint  exhibits  flexural  deformation  under 
tensile  loading  due  to  the  misalignment  between  the  neutral 
surface  and  the  loading  axis.  This  eccentricity  induces  a 
moment  distribution  along  the  joint  (see  Figure  9)  which 
acts  to  align  the  neutral  surface  with  the  loading  axis. 

The  analysis  of  the  "joggle-lap”  joint,  shown 
previously  in  Figure  1,  is  an  extension  of  the  concept 
discussed  above.  Again  it  is  clear  that  under  tension 
the  joint  will  experience  a  lateral  deflection  as  the 
neutral  axis  attempts  to  align  with  the  applied  force.  To 
analyze  the  joint  behavior  under  tensile  loading  conditions, 
it  was  decided  to  divide  the  joint  into  five  segments.  The 
obvious  places  to  divide  the  joint  are  illustrated  in 
Figure  6  along  with  the  corresponding  identifying  labels  and 
global  coordinate  system.  Reference  to  beam  segments  via 
their  identifying  numbers  will  be  utilized  throughout  the 
remainder  of  this  analysis. 

In  general,  the  goal  of  the  analysis  will  be  to 
determine  the  displacements  of  the  neutral  axis  as  measured 
perpendicularly  from  the  undeformed  neutral  surface.  Once 
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FIGURE  6:  PIECEWISE  REPRESENTATION 
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the  deflections  are  known,  one  may  calculate  a  moment 
distribution  along  the  joint  and  thus,  determine  the  stress 
distribution  at  any  given  cross-section. 

The  initial  intention  of  such  an  investigation  was 
to  develop  a  closed  form  solution  for  the  stresses  within 
the  joint.  This  effort  was  soon  thwarted  by  the  non-linearities 
encountered  in  the  governing  equations  for  the  beam  elements. 
These  non-linearities  result  from  a  coupling  between  the 
moment  and  deflection  solutions,  as  will  be  evident  later. 

As  an  alternative  solution,  the  displacement  field  was 
obtained  via  numerical  integration  routines . 

Linear  elastic  beam  theory  states  that  for  a  beam 
under  general  loading  conditions ,  the  local  radius  of 
curvature  is  given  by 

r  -  EI  (3) 

R  ~  “M 

where  R  =  radius  of  curvature 

E  =  modulus  of  elasticity 

I  =  moment  of  inertia  about  the  neutral  surface 
M  =  applied  moment 

The  radius  of  curvature  may  be  written  in  terms  of 
the  lateral  deflection  as  given  by  Eq.  (4) 


d2y 

dx2 


u  +  (^)2i3/2 


(4) 
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Realizing  that  under  the  assumptions  made  with  regard 
*  dv  2 

to  small  deflection  theory,  the  term  (^)  will  be  negligible 
when  compared  to  unity.  Thus  one  arrives  at  the  governing 
equation  for  straight  beam  elements . 


d2y  _  M 
dx2  El 


(5) 


Since  the  material  system  is  relatively  stiff,  it  is  assumed 
that  small  deflection  beam  theory  will  yield  sufficiently 
accurate  results.  Thus,  one  may  write  a  governing  differential 
equation  for  each  segment  of  the  joint.  By  matching  boundary 
conditions  of  deflection  and  slope  at  each  interface,  the 
deflection  of  the  entire  joint  may  be  obtained  as  a  function 
of  distance  along  the  neutral  axis.  Details  of  the  analysis 
may  be  referenced  in  Appendix  B. 

To  enhance  one's  understanding  of  the  joint  behavior 
under  applied  tensile  loadings,  Figures  7  through  10  show 
deflection,  slope,  moment,  and  shear  diagrams  respectively 
at  a  load  of  200  lbs.  Many  of  the  discontinuities  found 
in  the  plots  arise  from  a  shift  in  the  neutral  axis  which  is 
a  common  occurrence  among  lap  joints. 

It  was  stated  previously  that  analyzing  the  "joggle- 
lap"  joint  under  tension  was  a  non-linear  problem.  This  was 
seen  by  the  fact  that  the  moment  was  a  function  of  the 
deflection.  Another  way  to  view  the  non-linearities  of  the 
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FIGURE  7 :  DEFLECTION  OF  NEUTRAL  SURFACE  AT  FAILURE  LOAD 
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FIGURE  8:  SLOPE  OF  THE  NEUTRAL  AXIS  AT  THE  FAILURE  LOAD 
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joint  behavior  is  to  investigate  the  response  of  the  joint 
to  varying  tensile  loads.  Figures  11  and  12  provide  a 
clear  indication  of  the  deviation  from  linearity  even  for 
small  values  of  load.  Both  the  deflection  (Figure  11)  and 
moment  (Figure  12)  were  recorded  at  the  beginning  of  SEG3 . 
(i.e.  S3  =  0) 


Moment  (in-lbs)  Deflection  (in) 
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RESPONSE  OF  THE  JOGGLE-LAP  JOINT  AT  Sj  =  0 


b.  Finite  Element  Model  (tension) 

Anticipating  the  shortcomings  of  a  beam  bending 
model  in  the  adhesive  zone,  defined  to  be  the  area  of  actual 
bonding,  it  was  decided  to  model  this  area  using  finite- 
element  methods.  One  of  the  underlying  assumptions  of  small 
deflection  beam  theory  is  that  plane  sections  remain  plane 
during  pure  bending  action.  Clearly  the  validity  of  this 
assumption  is  questionable  in  the  bonded  area.  Another 
reason  for  employing  the  finite  element  technique  was  to 
uncover  any  local  stress  concentrations  that  may  not  be 
revealed  in  a  beam  analysis.  The  finite-element  mesh, 
consisting  of  7  material  types,  is  shown  in  Figure  13. 

Boundary  conditions  in  the  form  of  concentrated  loads  were 
applied  to  each  of  the  finely  meshed  ends.  Loading  conditions 
were  applied  away  from  the  adhesive  layer  at  a  distance  of 
1.5  times  the  thickness  in  an  effort  to  minimize  the  effects 
of  the  end  loads  upon  the  stress  solution.  An  explanation 
of  how  these  boundary  conditions  were  determined  will  follow 
shortly.  A  plane  stress  analysis  was  utilized  to  calculate 
the  displacement  and  stress  fields.  Figures  15  through  17 
are  the  result  of  a  plotting  routine  which  displays  lines 
of  constant  stress.  The  figures  should  be  interpreted  in  the 
same  manner  as  that  of  a  topographical  map.  Adjacent  lines 
spaced  closely  together  indicate  areas  of  high  stress 
gradients  and  possible  sites  for  structural  failure.  The 
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FIGURE  15:  THE  ADHESIVE  ZONE  IN  TENSION  -  CONTOURS  OF  a,  STRESS 


FIGURE  16:  THE  ADHESIVE  ZONE  IN  TENSION  -  CONTOURS  OF  a~  STRESS 


FIGURE  17:  THE  ADHESIVE  ZONE  IN  TENSION  -  CONTOURS  OF  t,„  STRESS 
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figures  are  labeled  according  to  the  component  of  stress 
being  displayed.  All  three  plots  are  the  result  of  loading 
the  specimen  at  the  tensile  failure  load  and  are  representative 
of  the  deformed  geometry. 

The  limitations  of  the  beam  bending  model  are 
clearly  displayed  in  Figure  15  and  reveal  the  justification 
for  the  finite-element  model.  Shown  in  the  figure  is  a 
smooth  transition  of  stress  across  a  change  in  cross-sectional 
area,  (i.e.  shift  of  the  neutral  axis)  as  calculated  by  the 
finite-element  method.  Experimental  results  have  shown 
this  to  be  a 'correct  representation  of  the  stresses.  Beam 
analysis  would  have  shown  a  sharp  discontinuity  in  the  stress 
profile  where  such  a  shift  in  the  neutral  axis  occurs.  Since 
the  moment  is  nearly  constant  throughout  SEG4  (see  Figure  9) 
beam  analysis  would  calculate  stress  contours  parallel 
to  the  adhesive  layer.  The  a a2,  and  t^2  stress  components 
are  global  oriented  stresses  as  opposed  to  those  that  can  vary 
according  to  element  orientation.  Marked  on  each  figure 
are  those  areas  where  the  assumptions  made  via  beam  analysis 
quite  appreciably  affect  the  accuracy  of  a  correct  solution. 

Many  analyses  of  lap  joints  assume  a  condition  of 
constant  shear  stress  in  the  adhesive  layer  itself.  This 
would  indeed  be  the  case  if  the  adherents  were  infinitely 
stiff  as  compared  to  the  adhesive  and  also  if  the  existence 
of  a  load  transfer  area  was  prohibited.  Shear  stress  data 
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from  the  finite  element  model  Is  plotted  in  Figures  18  and  19 
and  the  indication  is  clear  that  the  shear  stress  is  not  a 
constant  in  the  load  transfer  area.  The  case  of  constant 
shear  stress  found  toward  the  center  of  the  adhesive  zone, 
however,  reveals  the  linear  nature  of  the  displacement 
function  through  the  adhesive  thickness  in  this  area. 
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FIGURE  18:  SHEAR  STRESS  VARIATION  THROUGH  THE  THICKNESS  OF  THE 
ADHESIVE  ZONE 


FIGURE  19: 


Shear  Stress  (psi) 


SHEAR  STRESS  VARIATION  THROUGH  THE  THICKNESS  OF  THE 
ADHESIVE' ZONE 
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Boundary  Conditions  for  the -Finite-Element  .Model 

The  boundary  conditions  for  the  finite-element  model 

are  determined  by  applying  the  stress  distribution  as  directed 

by  the  beam  bending  model  to  the  finely  meshed  ends  of  the 

undeformed  geometry  of  the  finite  element  model.  In  other 

words,  the  stresses  in  the  deformed  geometry  (beam  model) 

must  be  moved  through  a  distance  to  their  equivalent  point 

of  application  in  the  undeformed  geometry  (finite-element 

model) .  The  reason  for  this  difficulty  with  boundary 

conditions  is  that  we  are  currently  utilizing  a  linearized 

■  2 

finite  element  routine,  SAP  V  ,  to  solve  a  non-linear  problem. 
Justification  of  such  a  procedure  will  hopefully  become 
lucid  with  time. 

To  facilitate  the  derivation  of  a  transformation 
routine,  Figures  20  and  21  illustrate  the  following  sign 
conventions.  Figure  20  depicts  a  stress  distribution  for 
the  left  hand  face  of  the  finite  element  model  with 
tension  being  taken  as  positive  and  compression  being 
negative.  Note  that  the  neutral  axis  is  not  coincident 
with  the  centroidal  axis  inherent  in  the  analysis  of  a 
curved  beam.  As  mentioned  previously,  this  fact  yields  a 
hyperbolic  stress  distribution  which  slightly  complicates 
the  computations.  (SEE  derivations  of  governing  equation 
for  stresses  in  a  curved  beam.  Appendix  A) 


2 

Structural  Analysis  Program  V;  University  of  Southern 
California,  Department  of  Civil  Engineering,  Oct.  77. 
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Figure  21  reveals  a  planar  yiew  of  the  deformed 
and  undeformed  sections.  It  is  assumed  in  this  derivation 
that  the  section  of  the  beam  can  at  most  undergo  a  trans¬ 
lation  and  a  rotation.  Translations  are  measured  via  the 
parameter  DEFLA  and  are  positive  radially  outward  as  shown. 
Small  deflection  theory  also  allows  the  rotations  to  be 
written  as  a  change  in  slope.  This  parameter  is  DUDSA  and 
is  positive  counter-clockwise. 

With  these  sign  conventions  clearly  in  mind  the 
stress  distribution  of  the  deformed  geometry  may  now  be 
resolved  into  concentrated  force  components.  Representing 
the  hyperbolic  stress  distribution  as  equivalent  point 
forces  and  point  couples  acting  at  nodal  points  labeled 
1  through  5  on  Figure  21  corresponds  mathematically  to  an 
integration  of  the  stress  distribution  between  fixed  limits. 


m 


M 

u  a 


/hi 
hi-l 


u  du  +  j. h±  F  cos  (e  +  DUDSA)  du  (6) 


R-u 


hi-l 


where  i  =  1-5 

Fn^  =  nodal  force  component 
M  =  moment 

u  =  distance  between  neutral  and  centroidal  axes 
a  =  cross-sectional  area 
R  =  radius  of  curvature 
F  =  load 

0  =  angle  subtended  by  SEG3 
DUDSA  =  local  slope  of  deformed  neutral  axis 
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The  first  term  of  Eq.  (.6).  represents  the  contribution  from 
the  hyperbolic  stress  distribution.  The  second  term  acts 
to  superimpose  the  component  of  force  due  to  longitudinal 
loading. 

A  correcting  moment  is  calculated  for  each  node  to 
equilibrate  the  two  representations  of  stress  on  the  section. 


M 


corr 


a (u) u  du  -  Fn^u 


(7) 


The  need  for  the  correcting  moment  is  due  to  the  fact  that 
a  distributed  force  is  now  represented  by  a  point  force  as 
shown  in  Figure  22. 

The  next  step  follows  from  a  translation  of  the 
point  forces.  Elementary  statics  dictates  that  a  point 
force  may  be  equivalently  represented  by  the  same  point 
force  and  an  added  moment  to  account  for  the  translation 
from  the  original  line  of  action. 

After  carrying  out  a  similar  procedure  for  the 
stresses  at  the  right  hand  side  of  the  finite  element  model, 
the  entire  system  is  set  in  equilibrium  by  accounting  for 
the  shear  acting  on  each  face  of  the  model.  The  values  of 
shear  are  obtained  directly  from  the  beam  bending  model. 

Thus  a  correct  set  of  boundary  conditions  has  been  determined 
for  the  finite-element  model  of  the  adhesive  zone.  A 
computer  routine  designated  by  CONVERT  was  written  to 
calculate  appropriate  boundary  conditions  and  may  be  found 


Methods  of  Analysis. 

B.  Flexure  Loading 
a.  Beam  Model 

The  bending  behavior  of  the  "joggle- lap"  joint  was 
also  studied.  It  was  found  that  the  theoretical  analysis 
was  far  simpler  than  that  encountered  for  tensile  loading. 

Each  segment  of  the  joint  (see  Figure  6)  was  modeled  as  if 
it  were  in  pure  bending.  Stresses  in  the  straight  beam 
numbers  were  calculated  via  the  flexure  formula  while  for  the 
curved  beams  the  formula 

My  _  (8) 

y  (R-y)ya 

where  M  =  moment 

y  =  coordinate  from  the  neutral  surface 
(positive  radially  inward) 

R  =  radius  of  curvature 

y  =  distance  between  centroidal  and  neutral  axes 
a  =  crossectional  area 

was  used. 

In  order  to  compute  bending  stresses  in  SEG4 
(layered  beam)  it  is  necessary  to  introduce  the  notion  of 
equivalent  sections.  In  this  method  we  assume  all  materials 
to  have  the  same  modulus  of  elasticity.  By  replacing  the 
actual  section  with  a  mechanically  equivalent  one  allows 
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the  flexure  formula  to  be  used  as  a  means  of  computing 
stresses.  The  width  of  the  sections  are  varied  so  that 
the  new  width  equals  the  ratio  of  the  old  modulus  of  the 
material  to  the  new  modulus  of  the  material  times  the  old 
width  as  shown  in  Figure  23.  Computing  I  for  the  specimen 

CvJ  • 

geometry, 


Ieq  =  I  ( 
i=l 


TZ  bihi 


3  +  ^j2) 


b  =  length  of  base 
h  =  length  of  side 
a  =  area 

d  =  distance  between  element  neutral  axis  and 
overall  section  neutral  axis 

it  is  apparent  that  the  effect  of  the  adhesive  layer  on 

overall  section  stiffness  is  negligible.  Using  the  flexure 

formula  and  the  relation 


(ax}  actual 


Jold 


new 


(a  ) 

'  x'equiv. 


the  stresses  in  SEG4  may  easily  be  calculated. 


b.  Finite  Element  Model  (flexure! 

The  boundary  conditions  of  the  finite-element  model 
may  be  changed  to  accommodate  pure  bending.  By  utilizing 
couples  at  the  finely  meshed  ends  of  the  model,  stresses  in 
the  adhesive  zone  may  be  monitored  where  it  has  been  shown 
that  the  results  from  beam  theory  are  less  accurate.  Figures 


24  through  26  display 


and  t^2  stress  contours 


respectively  within  the  "joggle-lap"  joint  in  pure 


bending . 
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FIGURE  24:  THE  ADHESIVE  ZONE  IN  BENDING  -  CONTOURS  OF  o,  STRESS 


FIGURE  25:  THE  ADHESIVE  ZONE  IN  BENDING  -  CONTOURS  OF  o0  STRESS 


FIGURE  26:  THE  ADHESIVE  ZONE  IN  BENDING  -  CONTOURS  OF  STRESS 


IV.  Experimental  Results 
A.  Tension 

As  set  forth  in  the  objectives  of  such  a  study, 
an  emphasis  was  to  be  placed  upon  developing  joint 
geometries  which  will  accommodate  high  rate  fabrication 
techniques.  In  an  effort  to  meet  this  criterion 
experimentally,  it  was  necessary  to  utilize  a  joint 
configuration  currently  being  molded  in  industry.  The 
time  and  expense  of  developing  in-house  molding 
capabilities  proved  to  be  beyond  the  scope  of  the 
research  at  hand.  Thus,  test  sections  were  cut  from 
premolded  panels  of  SMC  which  were  later  bonded  together 
to  form  the  joint. 

The  bonding  operation  was  also  directed  toward 
high  fabrication  procedures.  All  test  specimens  were 
adhesively  joined  at  Goodyear  Adhesives  Division,  Ashland, 
Ohio,  via  production  adhesives  application  techniques. 

It  was  felt  that  by  using  these  sophisticated  application 
procedures  optimum  adhesive  properties  could  be  obtained. 

In  general,  SMC  is  defined  to  be  an  anisotropic 
material  because  of  the  substantial  difference  between 
in-plane  and  out-of-plane  properties.  Referring  to  the 
coordinate  system  of  Figure  1,  the  constitutive  relations 
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FIGURE  27:  ELASTIC  RESPONSE  DUE  TO  TENSION  -  TOP  FIBER  STRESSES 
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ELASTIC  RESPONSE  DUE  TO  TENSION  -  BOTTOM  FIBER  STRESSES 


that  the  joint  invariably  strained  beyond  the  small-deflection 
range  at  considerably  small  loadings.  It  was  therefore  a 
rather  arduous  task  to  approximately  determine  the 
experimentally  applied  moment  to  the  joint.  The  correlation 
between  the  theoretical  and  experimental  data  may  be 
referenced  in  Figures  31-34.  As  in  the  case  of  tensile 
loading,  it  should  be  noted  that  the  stresses  in  SEG1  are 
again  considerably  higher  than  those  predicted  by  theory, 
which  is  attributable  to  the  molded  geometry. 


FIGURE  31:  TOP  FIBER  STRESSES  DUE  TO  BENDING 
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FIGURE  33:  BOTTOM  FIBER  STRESSES  DUE  TO  BENDING 
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FIGURE  34:  BOTTOM  FIBER  STRESSES  DUE  TO  BENDING 


Experimental  Results  from  Tension 
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*after  failure  initiation,  it  was  observed  that  the  crack  was 
propagated  via  interlaminar  shear 


V.  Failure  Analysis 

One  of  the  most  important  parameters  to  predict 
in  a  study  of  this  type  is  the  ultimate  loading  conditions. 

This  in  essence  dictates  the  choice  of  a  failure  criterion. 

The  maximum  stress  theory  will  be  employed  in  this  report 
because  of  its  simplicity  in  application  and  execution. 

Other  popular  failure  criteria,  such  as  the  Tsai-Wu 
criterion  were  deemed  inappropriate  due  to  the  limiting 
assumptions  made  in  accordance  with  beam  theory. 

Maximum  stress  criterion  states  that  the  material 
will  fail  when  any  component  of  stress  exceeds  the 
corresponding  material  strength.  In  general,  the  above 
statement  may  be  written  in  equation  form  as 

(14) 

(15) 

(16) 

T 

where  =  ultimate  tensile  strength 

C 

=  ultimate  compressive  strength 
=  ultimate  shear  strength 

These  equations  simplify  to  those  listed  below  after 
employing  the  local  coordinate  nomenclature  for  the 
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" joggle-lap"  joint. 
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App lying  this  failure  criterion  to  the  model,  it  was 
found  that  the  bottom  fiber  tensile  stresses  (see  Figure  35) 
predicted  the  ultimate  loading  of  the  joint  within 
experimental  error.  Thus  the  maximum  flexural  stress 
was  utilized  to  predict  failure. 

All  failures  occurring  as  a  result  of  tensile 
loading  were  initiated  along  the  bottom  surface  of  SEG3 . 

Crack  initiation  was  observed  to  be  of  the  net  tension  mode, 
while  propagation  appeared  to  be  due  to  "interlaminar  shear". 
There  was  a  general  consistency  among  the  initiation  and 
propagation  of  the  crack  for  all  tension  tests. 

It  was  thought  at  one  time  that  the  curved  sections 
of  the  joint  (SEG2,  SEG3)  were  either  fiber  deficient 
or  highly  anisotropic  yielding  a  potential  low  strength 
area.  However,  a  photomicrograph  of  this  cross-sectional 
area  clearly  shows  no  such  tendencies.  (See  Plate  7) 


(ksi) 
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FIGURE  35:  BOTTOM  FIBER  STRESSES  AT  THE  FAILURE  LOAD 


(ksi) 
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Ultimate  Tensile  Stresses 


Top  Fiber  Stresses 

_  Theoretical 

•  Experimental 


Ultimate  Compressive  Stresses 


FIGURE  36:  TOP  FIBER  STRESSES  AT  THE  FAILURE  LOAD 


VI.  Conclusions 


The  response  of  the  "joggle-lap"  joint  was 
investigated  for  both  tensile  and  bending  loads  in  this 
report.  It  was  found  that  experimental  data  correlated 
rather  well  to  the  values  of  stress  predicted  by  the 
analytical  model.  The  results  of  the  bending  study  were 
not  as  favorable ,  in  that  experimental  verification  proved 
to  be  more  difficult. 

A  parametric  study  was  undertaken  for  the  "joggle- 
lap"  joint  subject  to  tensile  loads  in  an  effort  to 
isolate  the  crucial  design  parameters.  In  Figures  37 
through  40  a  normalized  stress  value  is  plotted  against 
one  of  four  parameters  -  adherent  thickness,  inside  radius, 
contact  area,  and  load.  From  these  design  curves  the 
following  conclusions  are  inferred. 

•  If  weight  saving  requirements  are  not 
stringent,  the  effect  of  increasing 
adherent  thickness  drastically  reduces 
maximum  flexural  adherent  stress. 

•  Increasing  the  radius  of  curvature 
has  a  negligible  effect  on  reducing 
maximum  adherent  stress  due  to  a 
trade-off  between  mechanisms. 
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FIGURE  37:  EFFECTS  OF  ADHERENT  THICKNESS  ON  JOINT  STRENGTH 
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FIGURE  38:  EFFECTS  OF  INSIDE  RADIUS  OF  JOINT  STRENGTH 
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FIGURE  39:  EFFECTS  OF  CONTACT  AREA  ON  JOINT  STRENGTH 
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FIGURE  40:  EFFECT  OF  LOAD  ON  JOINT  STRENGTH 
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•  Neglecting  local  stress  concentrations, 
the  effect  of  reducing  the  overlap 
length  does  not  increase  adherent  stress 
significantly . 

•  In  the  -region  of  the  failure  load,  the 
maximum  adherent  stress  increases  linearly 
with  load. 

An  important  parameter  in  joint  design  is  that  of 
joint  efficiency.  This  parameter  is  defined  to  be  the 

p 

ratio  of  ultimate  joint  load  divided  by  the  ultimate  load 
carried  by  the  material  if  the  joint  were  not  present. 

The  joint  efficiency  of  the  "joggle-lap"  joint  in  tension 
is  calculated  to  be  0.153. 

The  adhesive  system  employed  in  this  report 
proved  to  be  quite  adequate  from  a  structural  point  of 
view.  For  the  given  overlap  length  of  1  in  (2.54  cm) 
there  were  no  recorded  failures  in  the  adhesive  layer. 
Failure  loads  were  predicted  using  the  maximum  flexural 
stress  as  the  limiting  criterion. 

This  report  would  be  incomplete  if  it  did  not  offer 
several  suggestions  for  future  work  as  an  outgrowth  of 
this  study.  An  obvious  limitation  to  the  work  reported 
herein  is  the  inability  to  extensively  verify  the  analytical 
model  by  experimental  testing  of  various  joint  geometries. 
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Further  development  in  this  area  would  greatly  increase 
the  reliability  of  the  computer  model. 

More  detailed  work  needs  to  be  completed  in  the 
response  of  the  "joggle-lap"  joint  to  bending  loads. 

This  report  included  only  a  cursory  investigation  of 
bending  behavior  as  a  means  of  identifying  the  underlying 
problems  associated  with  the  experimental  verification 
of  theory. 

It  is  felt  that  this  report  will  provide  a 

fundamental  basis  for  future  research  concerning  the 

0 

"joggle-lap"  joint. 
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X.  Appendices 
Appendix  A 


Derivation  of  the  Governing  Equations 
for  a  Curved  Beam 


Consider  the  curved  beam  element  shown  in 


Figure  41. 


FIGURE  41:  CURVED  BEAM  ELEMENT 


The  analysis  begins  by  seeking  an  expression  for  the 
strain  distribution  perpendicular  to  the  neutral  axis. 
Assume  that  the  curved  beam,  with  an  initial  radius  of 
curvature  R,  undergoes  a  small  elastic  deformation  due 
to  the  applied  moment,  (.It  is  important  to  note  that 
the  neutral  axis  of  bending  for  a  curved  beam  does  not 
necessarily  coincide  with  the  centroidal  axis  of  the 
beam.)  Under  the  action  of  this  moment  it  becomes 
apparent  that  segment  cd  rotates  about  the  neutral  axis 
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to  a  new  position  c'd'.  It  is  assumed  here,  as  in 
classical  beam  analysis,  that  plane  sections  remain 
plane.  It  is  readily  seen  that  while  the  deformation 
of  the  beam  varies  linearly  with  the  distance  from  the 
neutral  axis ,  the  strains  do  not .  The  reason  is  that 
the  original  length  of  all  the  fibers  prior  to  the 
application  of  the  moment  are  not  constant. 

Thus  the  following  relation  for  the  strain 
distribution  is  written  below. 


e£  ^  -yA0 
A7  (R-y)  A<j> 


(20) 


where  e^  =  elongation 

y  =  radial  coordinate  (positive  radially  inward) 
A0  =  angle  of  deformation 
Acj>  =  angle  subtended  by  curved  beam 
The  above  equation  shows  the  strain  to  vary  hyperbolically 
across  the  section.  Using  the  plane  stress  constitutive 
relation,  Eq.  (20)  becomes 


„  _  -Ey  A0 
(R-y)  A(() 


(21) 


Now  it  is  appropriate  to  derive  the  formulas  for  flexural 
stress.  First  assume  that  the  portion  of  the  beam  is  in 
equilibrium.  Following  directly  one  may  write  the 
equations  of  equilibrium  for  an  arbitrary  section. 
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IF  .  .  = 

axial 

/Aada  =  0 


0 


Making  the  appropriate  substitutions  for  the  stress 
Eq.  (22)  becomes 


(22) 


-EyA6da 
A  (R-y)A<j> 


(23) 


Assuming  E,  A<j>,  and  A0  to  be  constants  the  integral  is 
simplified  as  shown  in  Eq.  (24) . 


r  yda 

ja  Tr^yT 


0 


(24) 


It  is  possible  to  solve  Eq.  (24)  for  the  radius  of  curvature 
and  thus  locate  the  neutral  surface;  however,  it  will  suffice 
to  let  Eq.  (24)  stand  as  is  for  now.. 

Referring  to  Figure  42  and  summing  moments  about  the 


section 


neutral  axis,  one  finds  that  the  stress  distribution  must 
also  satisfy  the  equation  below. 
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M  =  -/A a yda 


(25) 


making  the  appropriate  substitutions,  Eq.  (25)  becomes 


M  -  -'a 


(26) 


„  _  EA0  r  y 
M  Act)  fA  7r Iy)d 


(27) 


Notice  the  algebraic  relation  that  permits  the  substitution 
of  an  equivalent  expression  into  Eq.  (27) . 


Y  =  P-y 

R-y  R-y 


-  y 


(28) 


Eq.  (27)  now  becomes 


EA9  rj-  P yda - j.  yda] 

w  Ac(> L  J  A  R-y  A  *  J 


(29) 


and  from  the  result  of  Eq.  (24) 


M  =  (  R ( 0 )  -  au) 


(30) 


where 


a  =  area 

u  =  distance  between  the  neutral  and 
centroidal  axes 


Rearranging  Eq.  (30)  yields 


A0 

A<t> 


-M 

Eau 


(31) 
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Comparing  this  equation  with  the  well-known  deflection 
equation  for  straight  beams,  it  is  apparent  that 

dfz  =  Ji  (32) 

dx2  El 

the  left  hand  side  of  Eq.  (31)  is  not  yet  suitable.  The 
ultimate  goal  of  such  an  analysis  is  to  seek  an  equation 
that  relates  the  deflection  of  the  neutral  axis  to  the 
position  along  the  neutral  axis. 

Consider  Figure  43  shown  below. 


The  beam  is  deflected  as  shown  to  illustrate  the  most 
general  case  of  a  non-constant  moment.  That  is,  the  moment 
is  a  function  of  position.  Now  the  deflection  can  be  measured 
as  the  deviation  between  the  undeformed  neutral  surface  and 
the  deformed  neutral  surface.  For  convenience  just  the 
neutral  axis  and  appropriate  parameters  are  drawn  in 
Figure  44. 
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A  coordinate  system  u,s  is  defined  and  shown  in  the  figure 
where  s  traverses  tangentially  to  the  undeformed  neutral 
axis  and  u  is  defined  to  be  perpendicular  to  that  axis. 

Enlarging  the  area  of  interest  and  focusing  on 
the  triangle  of  Figure  45,  one  finds  that 
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9  (s) 


Au 

AA 


Realizing  that  tan  a  =  a  for  small  a,  it  follows  that 


AA 


(R+u) AS 
R 


and  thus 


6 (s)  =  RAu 

(R+u) As 


which  may  be  written  as 

Au  _  (R+u)  9  (s) 

As  R 

Finally  in  the  limit  as  As->-0 :  Eq.  (33)  becomes 

lim  Au  _  (R+u) 9 (s)  _  du 

As-*0  As  R  ds 

From  Eq.  (31)  ,  several  simplifications  can  be  made  with 
the  proper  substitutions . 


(33) 


(34) 


A0 

A0 


-M 

Eau 


where 


Acf> 


As 

R 


lim  A0 
As-*-  0  As 


-M 

REau 


d0 

ds 


(35) 


Differentiating  Eq.  (34)  with  respect  to  s  yields 
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d2u  _  (R+u)  d6 

ds^  R  ds 


(36) 


and  substituting  Eq.  (35)  into  Eq.  (36)  yields  the  final 
results  -  a  second  order  differential  equation  relating 
deflection  to  position  in  terms  of  the  applied  moment. 


d2u  _  ( R+u ) M 

2  o  - 

ds  R  Eau 


(37) 


Appendix  B 


Beam  Bending  Model  of  the  "Joggle-Lap"  Joint 
SEG1  may  be  modeled  as  a  straight  beam  shown  in 
Figure  46. 

.  y 


FIGURE  46:  'SEG1  MODELED  AS  A  STRAIGHT  BEAM 


In  general,  the  moment  experienced  by  any  segment  originates 
from  two  sources :  eccentricity  from  geometry  and 
eccentricity  due  to  deflection.  The  preceding  statement 
may  be  written  algebraically  as  follows. 


where 


M  =  F  (e 


geom 


edefl} 


M  =  moment 
F  =  applied  force 
e  =  eccentricity 


(38) 


It  is  readily  seen  that  egeom  =  0  for  SEG1.  Writing  Eq.  (38) 
in  the  local  coordinate  system,  the  moment  experienced  by 
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this  segment  reduces  to 


M  =  FUj^  (39) 

where 

=  deflection  in  the  local  coordinate  system 
Substituting  Eq.  (39)  into  Eq.  (5)  yields 

,2 

d  u,  Fu, 

— ~  — —  =  o  (40) 

ds^  El 


The  corresponding  boundary  conditions  are  expressed 
below 


where  u  is 

o 

The  solution 
to  be 


ux(0)  =  0 
ul(£l)=  uo 


yet  undetermined. 


of  Eq.  (40)  is  of  standard  form  and  known 


U1  =  C1  s^n51  V  F/EI  •  s]_  +  C2  cosh  v  F/EI  •  sx  (41) 

Applying  the  boundary  conditions  to  Eq.  (41)  determines 
the  constants  and  C7  to  be 


0 

uQ  /  sinh  / F/EI  • 


and  thus 
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u,  =  Uo  sinh  /F/EI  •  s. 

1  —7= - i  0<s, <£,  (42) 

sinh/  F/EI  •  A.  A  x 

where  u0  is  necessarily  negative  to  correspond  with  the 
physical  system.  In  other  words,  for  a  given  tensile 
load  it  is  expected  that  SEG1  will  deflect  downward. 

(Figure  6) .  Also 


_  _ 

(A,)  =  u0  /  P/EI  cosh  7  F/EI  A. 


sinh  /  F/EI  A. 


It  should  be  noted  that  the  deflection  as  given 
by  Eq.  (42)  is  not  known  explicitly  in  terms  of  the  given 
parameters.  U0  is  still  unknown  and  it  will  be  shown  later 
how  this  value  may  be  determined  uniquely. 

SEG2  is  modeled  as  a  curved  beam  and  shown  in 
Figure  47.  The  local  coordinate  system  is  a  curvilinear 
coordinate  system  with  the  s2  axis  traversing  the  neutral 
axis  as  shown.  Positive  deflections  are  measured  normal 
to  the  undeformed  neutral  axis  in  the  direction  of  u2 . 


centroidal  axis 
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From  the  derivation  of  the  general  case  for  a  curved  beam 
in  pure  bending  (see  Appendix  A) ,  the  governing  equation 
for  the  deflection  is 

d2u-  (R  +  u9)M 

- £  =  — _ — -±_  (43) 

ds^  R  Eau 

where 

S2  =  arc  length 

U2  =  deflection  normal  to  neutral  axis 
M  =  moment 

R  =  radius  of  curvature 

E  =  modulus  of  elasticity 

a  =  cross  sectional  area 

u  =  distance  between  neutral  axis  and 
centroidal  axis  and  its  value  is 
necessarily  negative 

The  moment  may  be  written  as  the  product  of  the  applied 
load  and  the  eccentricity,  where  the  eccentricity  in  this 
case  consists  of  both  geometry  and  deflection  considerations. 

At  this  point,  it  is  appropriate  to  introduce 
the  notion  of  extensional  effects.  It  is  realized  that 
with  the  given  loading  conditions,  the  "joggle- joint "  will 
undergo  deflections  parallel  to  the  neutral  axis  as  well. 

This  fact  would  be  of  little  concern  if  all  beam  segments 
of  the  joint  configuration  had  their  neutral  axis  aligned 
with  the  loading  axis.  If  this  were  the  case,  the 
longitudinal  displacement  would  not  affect  the  eccentricity. 
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However,  it  is  evident  that  the  extensional  strains  in 
the  curved  beam  segments  give  rise  to  an  added  component 
of  eccentricity  defined  to  be  e  To  calculate  the 

value  of  e  t,  one  merely  applies  the  criterion  of 
force  equilibrium  to  SEG2  (Figure  47)  in  the  local 
coordinate  system. 


IF  =0 
u2 

ZFa  =  0 
S2 


H2cos0  +  V2sin0  =  F 


H2sin0  =  V2cos0 


H2  =  Fcose 


thus 

where  0  =  angle  subtended  by  SEG2 

Employing  the  constitutive  relationship 
a  =  Ee 

where  a  =  stress 

E  =  modulus  of  elasticity 
£  =  strain 

and  considering  only  the  y  (global  coordinate)  component 
of  the  extension  we  thus  arrive  with  the  expression  for 

eext* 


eext  =  Fs2cos (s2/R) sin (s2/R) 


aE 


Eq.  (44)  must  be  added  to  the  other  terms  which  comprise 
the  eccentricity  due  to  deflection. 

Therefore  Eq.  (43)  becomes 


(44) 


R  Eau 


where 


e  _  =  R ( 1-cos  (-p)  +  u) 
geom  s  k 

edefl  "  u2cos  ‘“I* 


eext  =  Fs2cos  (s2/R)sin  (s2/R) 


Initial  conditions  for  SEG2  are  found  by  matching 
deflection  and  slope  at  the  1-2  interface. 


u2(0)  =  u0 


—A  (0)  =  u0  J  F/EI  cosh  sj  F/EI  Jl, 


sinh  v  F/EI 

Using  a  numerical  integration  routine  to  solve  Eq.  (45)  the 
deflection  u2  may  be  marched  out  as  a  function  of  arc  length 
s2.  A  Runge-Kutta  method  based  on  Verners  fifth  and  sixth 
order  pair  of  formulas  was  used.  An  explanation  of  the 
integration  routine  DVERK  may  be  referenced  in  Appendix  C. 

Figure  48  shows  SEG3  modeled  as  a  curved  beam. 

From  Eq.  (37)  the  governing  differential  equation  for  a 
curved  beam  in  pure  bending  is 


(r+u3) 

R2Eau 


where  M  =  F(egeom  +  edefl  +  eext) 
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egeom  =  ^  +  R(l-cos0)  -2ucos0  +  Rsin(^-  -9  +  s3/R)-cos9 


where  9  =  angle  subtended  by  SEG3 

u  =  distance  between  centroidal  and  neutral 
axes 

R  =  radius  of  curvature 

s,  =  arc  length  along  neutral  surface  of 
SEG3 

Also  edefl  =  M  cos  (ir/2  -  9  +  s3/R)  (48) 


From  a  similar  argument  developed  earlier  it  may 
be  shown  that  eext  for  SEG3  is  given  by 


eext  =  Fs3cos (0-s3/R)  sin  (0-s3/R)  (49) 

aE 


(47) 
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Thus  Eq.  (46)  becomes 


-(R+u)  F  (e  ^  +  e,  +  e  ,) 

~ —  '  geom  defl  ext 

R  Eay 


(50) 


where  egeQm,  edefl,  and  eQxt,  are  given  by  Eqs.  (47),  (48), 

and  (49)  respectively.  Matching  boundary  conditions  at  the 
2-3  interface  provides  initial  conditions  to  Eq.  (50)  which 
may  be  integrated  numerically  as  before. 

SEG4  is  analyzed  as  a  multi-layered  beam  and 
shown  in  Figure  49.  Treating  this  segment  to  be  composed 
of  three  linear  elastic  beam  elements,  the  governing 
differential  equation  follows  from  Eq.  (5)  with  a  slight 
modification. 


where 


M 


3  E  .  I . 
I  li 


i=l 


(51) 


1^  =  moment  of  intertia  of  the  ith  section 
about  the  neutral  axis 

E^  =  modulus  of  elasticity  of  the  ith  element 


3 

Z  E^I^  is  referred  to  as  an  effective  flexural  stiffness 
i=l 

and  is  merely  a  constant.  The  moment  is  defined  in  the 
usual  manner  as 

geom  +  edefl^ 


M  =  F(e 


FIGURE  49:  SEG4  MODELED  AS  A  LAYERED  BEAM 


FIGURE  50:  SEG5  MODELED  AS  A  STRAIGHT  BEAM 
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where 


egeom  =  -5<t+b) 
edefl  =  u4 


Thus  Eq.  (51)  becomes 


,2 

4  _  F  ( .  5  ( t+b)  +  U4) 

ds4  l  Vi 

i=l 


(52) 


Initial  conditions  are  found  by  equating  the  deflection 
and  slope  at  the  3-4  interface.  Following  in  the  usual 
manner,  Eq.  (52)  is  integrated  to  obtain  an  expression 
for  the  deflection  of  SEG4  as  a  function  of  arc  length 
in  the  local  coordinate  system. 

Finally  SEG5  is  shown  in  Figure  50  modeled  as 
a  straight  beam  member.  The  governing  differential 
equation  is  the  same  as  Eq.  (5) 


.2 

d  u. 


ds, 


M 

El 


(53) 


where  M  =  Fu^ 

and  the  initial  conditions  are  obtained  by  matching  the 
deflection  and  slope  at  the  4-5  interface.  Upon  integration 
of  Eq.  (53)  the  deflection  SEG5  will  be  a  known  function  of 
the  abscissa  s^  of  the  local  coordinate  system.  Therefore  the 
deflection  and  slope  at  point  P  of  Figure  50  are  also  known. 
But  it  should  be  apparent  that  the  values  of  the  deflection 


-93- 


and  slope  at  this  point  must  be  zero  or  at  least  within 

certain  tolerance  limits.  This  in  fact  is  the  final 

boundary  condition  to  the  problem  that  is  needed  to  uniquely 

determine  the  value  of  u0  which  was  previously  assumed  to 

be  arbitrary.  Thus,  through  an  iterative  process,  a 

correct  value  of  uQ  may  be  calculated  by  assuring  that  the 

3 

deflection  and  slope  of  point  P  of  SEG5  is  sufficiently 
close  to  zero.  To  avoid  confusion,  it  should  be  noted  that 
by  specifying  zero  deflection  at  point  P  we  will  force 
the  slope  to  zero  by  the  nature  of  the  deflection  function 
of  SEG5.  So  in  fact  this  is  a  well-posed  problem,  whereby 
we  specify  only  enough  boundary  conditions  as  there  are 
unknowns.  The  process  for  correctly  determining  uQ  is 
shown  schematically  in  Figure  51. 


3 

it  was  found  that  reliable  results  were  obtained  by 
using  the  tolerance  limits  listed  here. 

|deflection  (P) |  <  .00001 

| slope  (P) |  <  . 00005 


Deflection  Deflection  Deflection 
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Appendix  C 


Computer  Programs 
a .  JOGGLE 

To  facilitate  ease  of  calculation,  a  computer 
routine  identified  as  JOGGLE  was  developed  and  may  be 
referenced  below.  Essentially  this  program  calculates  a 
correct  value  of  u0  and  proceeds  to  determine  a  solution 
for  the  deflection  while  calculating  stress  profiles 
along  the  joint  configuration.  These  stress  profiles  are 
linear  in  the  straight  beam  members  (SEG1,  SEG5)  and 
hyperbolic  in  the  curved  beams  ( SEG2 ,  SEG3) . 
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C- 

U 

00000 

GGGGG  GGGGG  L 

EEEEE 

00000003 

c- 

J 

0  0 

G  G 

L 

E 

00000004 

c- 

J 

0  0 

G  GG  G 

GG  L 

EEE 

0  3000005 

c- 

J  u 

0  0 

G  G  G 

G  L 

E 

00000006 

c« 

UUJUU 

00000 

GGGGG  GGGGG  LLLLL 

EEEEE  ' 

00000007 

C- 

1 

01)000003 

c- 

0  3000009 

c- 

00000010 

c- 

0000001 1 

c- 

00000015 

c- 

0300001  6 

c- 

ANALYTICAL 

BEAM  BENDING  MODEL 

00000020 

c- 

FOR 

A  JOGGLE  LAP 

JOINT 

00000021 

c- 

0D000024 

c- 

03000025 

c- 

DEVELOPED  BY: 

RICHARD  C. 

GIVLER 

03000026 

c- 

UNIVERSITY  OF  DELAWARE 

C  3000027 

c- 

SEPT 

78  -  JAN  79 

03000028 

c- 

03000040 

c- 

03000041 

c- 

C  3000042 

c- 

C 3000043 

SSET 

AUTOBIND 

C 300 0100 

SRESET  FREE 

G  3000 1 10 

file 

6 ( K  IND  =  R  EMOTE . MAXRECSIZE= 

22) 

C  3000 1 20 

D I  MENS I  ON 

PROD ( 3  ) 

.  ERR  (  1  ) 

.  T ( 2 ) ,  HI ( 3  )  .  Y PRIME! 2 ) 

C 30001 50 

COMMON  R . 

PLOAD. 

ESMC,  WIDTH ,  THICK. 

YBAR.  THETA. 

EADH.  PI 

C 30001 75 

BOND.  J 

.  TORC. 

TORCS.  TRAC.  TRACE. 

SHEAR 

C 30001 76 

PARAMETERS  AND  NOMENCLATURE 


MATERIAL  THICKNESS 
TH I CK= . 1 

C- - LONGITUDINAL  MODULUS  OF  ELASTICITY  OF  SMC  IN  PSI 

£SMC=2 . 1 E+06 

C-- - MODULUS  OF  ELASTICITY  OF  ADHESIVE  IN  PSI 

E A  DH= 1  . OE+05 

C . LOAD  IN  LBS 

P  LOAD  =  20  0 . 

C . INSIDE  RADIUS  IN  INCHES 

RADI=2. 5* THICK 

C- - OUTSIDE  RADIUS  IN  INCHES 

RAD0  =  3 . 5  +  TH ICK 

C . BONOING  THICKNESS  IN  INCHES 

BOND= .03 

C----- CONTACT  WIDTH  IN  INCHES 
CONTA= 1 . 0 

C . SPECIMEN  WIDTH 

WIDTHS  .  0 

C- - LENGTH  OF  SEGMENT  1 

$EGA=3 . 5 

C- - - LENGTH  OF  SEGMENT  5 

SEGB=4 . 0 

C . LEFT  INTERVAL  LIMIT  FOR  ITERATION 

A I NT= - .04 

C-----RI GHT  INTERVAL  LIMIT  FOR  ITERATION 
B 1 NT=+ . 0 1 

C . TOLERANCE  LIMIT  ON  INITIAL  DISPLACEMENT 

ERR ( 1 ) = . 00001 


C 30001 78 
C3000180 
00000200 
C 00002  1C 
C 0000220 
COO 00300 
C  0000400 
C  0000500 
COOC06GG 
C  0000700 
C  0000800 
0 000090C 
COOOIOOO 
C0001 100 
C00012CG 
C  000 1 300 
C  0001 400 
CC001500 
C 000 1600 
C0001700 
COCO  1000 
C0001S00 
CC-002C00 
C  0002 1 00 
CC002200 
C0002300 
CC002400 
00002420 
C0002430 
CG002435 
00002440 
0C002445 
C0002450 
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c . STEP  SIZE  FOR  NUMERICAL  INTEGRATION 

STEP=50 . 

P I  =3 . 141592654 

C . TRACING  CONSTANTS 

TR4C= 1 . 0 
TR4CE=1 . 0 
C . 

C--  NOTE : DEF LECTIONS  ARE  MEASURED  NORMAL  TO  THE  UNDEFORMED 
C-  NEUTRAL  AXIS 

C- 

YO=( AINT+BINT )/2. 

C- 

c . ITERATIVE  CALCULATION  OF  YO  TO  FORCE  ZERO  DEFLECTION 

C . AND  SLOPE  AT  END  OF  SEG5 

C . yO  MUST  LIE  BETWEEN  THE  PROPOSED  LIMITS  OF  AINT  AND  BINT 

C- 

(J.****#*****^M**i******»**^^M***tM.»t#M**f*(*,T|^M¥ 

£.  ************  ************  ******************  *****  ■,****  ****** 

C  -  *  *  *  •-■******.  +  +  + 

C  NUMERICAL  INTEGRATION  VIA  LIBRARY  ROUTINES  DVERK  AND  UERTST 
C  •  *#*¥T^****J«!*if+  TV  »•****  *  at  A**********,,*****.,.**** 

£.**>  V  *#***************.*. 4.  *****.,, 

C-****v  *********************  .j,*,^*******  ****,,*****  ^^^^^^ 

c- 

101  CONTINUE 

DO  200  X=0 , SEGA+ .005 . SEGA/STEP 

PROD  (  3  }  =  SORT  (  PLOAD/  (  ESMC  M  W I  DTH'  TH I CK -  *  3  .  /  1  2  .  )  )  } 

PROD( 1  )  =  YO* (  . 5* ( EXPf PROD( 3 ) *X ) - EXP{ - PROD( 3 ) -  X ) ) ) 

PROD { 2 ) = . 5*(EXP(PR0D(3)"SEGA)-EXP{  -  PROD ( 3 ) *  SEGA )  ) 

T ( 1 ) =PROD ( 1 )/PR0D(2 ) 

T ( 2 ) =  Y  0 ' PROD ( 3 } /PROD ( 2 ) M  . 5*(EXP( PROD ( 3 ) *X )  +  EXP 
-f -PR0D(3)*X)  )  ) 

200  CONTINUE 

DIMENSION  C ( 24  )  ,  Y(2),  W{2.9) 

EXTERNAL  FCN1 

C- - CALCULATION  OF  THE  RADIUS  OF  CURVATURE  FOR  CURVED  MEMBERS 

R  =  WIDTH  *■  THICK  /  (  ALOGf  R  ADO/ RAD  I  )  ) 

NW  =  2 

C . CALCULATION  OF  THE  ANGLE  SUBTENDED  BY  SEG2  AND  SEG3 

TH£TA  =  ARCOS ( (5. *THICK -  BOND ) / ( 6 . -THICK)  ) 

N  =  2 

C . CALCULATION  OF  THE  DISTANCE  BETWEEN  NEUTRAL  AXIS  AND 

C- - CENTRO  I  DA  L  AXIS  OF  CURVED  MEMBERS 

Y3AR=RAD I+THICK/2. -R 
X  =  0.0 

Y  (  1  )  =  -  T  (  1  ) 

Y ( 2  )*-T( 2  ) 

TO  L=  . 00000 1 
IND-1 

DO  300  2=0.0. R*THETA+. 001 .R*THETA/STEP 
XEND=  FLOAT ( Z ) 

CALL  DVERK (N . FCN1 , X . Y , XEND . TOL , I ND . C . NW , W , I ER ) 

IF( IND. LT.O.OR. IER.GT.O)  GO  TO  20 
300  CONTINUE 
20  CONTINUE 
R I  NT  =  RNE  W 
EXTERNAL  FCN2 
X  =  0.0 
Y(  1  )  =  - Y (  1  ) 

Y  {  2  )  =  -  Y  (  2  ) 

NW  =  2 

N=  2 


0  300245 1 

03002452 

03002475 

00002500 

0000250 1 

00002502 

00002600 

00002700 

0  3002800 

03002900 

00003000 

00003100 

03003200 

00003210 

00003220 

0  3003400 

03002405 

00003410 

00003420 

C  0  0  0  3  4  2  5 

03003430 

00003435 

00003440 

00003445 

03003450 

00004200 

0000 4300 

00004400 

00004500 

00004600 

C  0004700 

00004800 

03005100 

00006000 

C  0006 1 00 

C0006150 

C 3006200 

C 0006300 

C0006350 

C0006400 

00006500 

C0006550 

C 3006552 

03006600 

C  0006700 

C  0006800 

C0006900 

C 3007000 

C0007100 

C0007200 

C  0007300 

C 0007400 

C0007500 

C 3007800 

C0007900 

C  0007950 

C  0008600 

C 0003700 

C  0008720 

C 0008730 

C00C8800 

C0008900 


o  o  o  o  o  o  o 
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IND=1 

DO  250  Ms  1  .24 
C I M ) =0 , 0 
250  CONTINUE 

DO  350  Z=0.0/R*THETA+.Q01 . R*THE7A/STEP 
XENO-FLO AT ( Z ) 

CALL  DVE  RK { N , FCN2 . X . Y , XEND , TOL. I ND . C , NW . W , I ER ) 

IFf IND. LT.O.OR. IER.GT.O)  GO  TO  70 
350  CONTINUE 
70  CONTINUE 

EXTERNAL  FCN3 
X  =  0.0 
NW  =  2 
N  =  2 
I  ND  =  1 

DO  290  Mai  ,24 
C(M)*0.0 

290  CONTINUE 

DO  400  Z=0.G.C0NTA+.005,C0NTA/STEP 
XEND=FLOAT(Z) 

CALL  DVE  RK ( N , FCN3.X . Y . XEND , TOL . I ND . C . NW . W , I ER ) 

IFf IND. LT.O.OR. IER.GT.O)  GO  TO  80 
400  CONTINUE 

80  CONTINUE 
EXTERNAL  FCN4 
X  =  0 . 0 

NW  =  2 
N  =  2 
I  ND*  1 

DO  291  M»1  ,24 
C  (  M  )  a  o  .  0 

291  CONTINUE 

00  246  Z=0 .0 . SEGB+ . 005 . SEGB/STEP 
X£ND=  FLO AT ( Z ) 

CALL  Dv'ERK  (  N  ,  FCN4  .  X  .  Y  .  XEND  .  TOL  .  I  ND  .  C  .  NW  ,  W  .  I  ER  ) 

IFf IND. LT.O.OR. IER.GT.O)  GO  TO  31 
246  CONTINUE 

81  CONTINUE 

IFf  ABSf  Y  (  1  )  )  .  LT.ERRM  )  )  GO  TO  88 
I F ( Y { 1 ) . GT . 0 )  GO  TO  100 
IFf Y ( 1 } . LT.O)  GO  TO  89 
100  CONTINUE 
81 NT=YO 

Y0=( AINT+BINT )/2. 

GO  TO  101 
89  CONTINUE 
A  I NT  =  YO 

Y0=( AINT+8INT )/2 . 

GO  TO  101 
88  CONTINUE 

_+*••«■**»•****»,******  ******  ****************** '*****,,  **£*.>£*** 
-  *  *  *  *  *  »»**»**  ******  ***************.■**************>  *  f  X  r  .  »  *  *  * 
„************  ********  *************  ********************** 


CALCULATION  OF  DEFLECTION  AND  SLOPE  AS  A  FUNCTION 
OF  X  FOR  SEGMENT  1 

WR I TE  f 6 . 500 )  YC 

500  FORMAT (//////////, 1  OX , r Y0= 1  , F 1 4 . 1 1  .  ) 

WR I TE ( 6 , 501 ) 

501  FORMAT  {  ' . - . '  ) 

WR I TE ( 6 . 502 ) 


0  30O9COO 
03009100 
00009200 
03009300 
03009400 
00009500 
03009600 
03009700 
0301 0000 
00010100 
0301 1000 
0301 1 100 
0001 1200 
0001 1 300 
0301 1400 
0301 1500 
0301 16C0 
03011 700 
0001 1800 
0301 1900 
03012000 
00012100 
0301 2400 
0301 2500 
030 1 3400 
03013500 
00013600 
0001 3700 
00013800 
C 301 3900 
03014000 
00014100 
03014200 
03014300 
03014400 
C301 4500 
C 301 4800 
C3014S00 
03014951 
COO  14962 
03014963 
0301 5C00 
03015075 
C  30 1  5076 
C3C1 5077 
03015080 
C 301 5085 
03015086 
COO  1 5087 
0001 5C39 
0001 6C20 
COO  16030 
CO0 16040 
C  3024000 
C002410G 
COO 2 4200 
C  0024300 
C  3024400 
C0024500 
C 3024600 
C  00247  00 
00024800 


no  o  o  o  o 
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502  FORMAT  (  7 . '  ) 

WRITE(6. 2) 

2  FORMAT {/////,  7  DEFLECTION  AND  SLOPE  FOR  SEGMENT  1  ',//)■ 

WRHEio.  9) 

9  FORMAT (  10X.  7X  DISTANCE '  .  5X  .  7 DEFLECTION  7  ,7X,  'SLOPE7  ,9X,  'STRESS' 
-.8X, 'MOMENT' ,//) . 
d  W  =  0 

DO  505  X  =  0 . 0 . SEGA* .005. SEGA/STEP 

PROD { 3 ) =  SORT ( P LOAD/ ( ESMC ' ! W I DTH ' THICK  *  * 3 . / 1 2 .  )  )  ) 

P ROD ( 1  )  =  YOM  .5* (EXP!  PROD ( 3 ) *X ) -  EXP { -3R0D(3)'X)  )  ) 

PROD (  2  )  = . 5* (EXP( PROD! 3) "SEGA) -EXP! -  PROD ( 3 ) *  SEGA ) ) 

T (  1  ) =  PRO  D ( t ) /PROD { 2 ) 

T ! 2  )  = YO ' PROD{ 3 )/PROD( 2  )*( . 5* ( EXP( PROD ( 3 ) *X )  +  EXP 
- (  -  PROD ( 3  )*X  )  )  ) 

C- - CALCULATION  OF  TOP  AND  BOTTOM  FIBER  STRESSES 

SIGT=-PL0AD*T(  1  )  *TH  I CK/2 .  /  {  1  . / 1 2 . *  W 1DTH* THI CK*  *3 )  +PLOAD/ 

- ( WIDTH  *  THICK ) 

S 1 GB=PLO AD*T ( 1  )*THICK/2./f 1  ./I  2. * WIDTH* THICK* *3 J  +  PLOAD/ 

- (WIDTH-THICK) 
d  =  dW+2 

IF(  INT(.J/2)  .NE.J/2.0)  GO  TO  199 
WRITE(6,10)  X.  T ( 1 ) .  T ( 2 ) .  SIGT.  PL0AD*T{1) 

10.  FORMAT (  T0X,5(E11  . 4 . 3X  )  ) 

WR I TE (  6  ,  11  )  SIGB 
.11  FORMAT { 5  2X  ,  E 1 1  .4.//) 

199  JWsJW+1 
505  CONTINUE 

----CALCULATION  OF  DEFLECTION  AND  SLOPE  AS  A  FUNCTION 
OF  ARC  LENGTH  F^R  SEGMENT  2 

WRITE! 6 . 24) 

24  FORMAT (/////, 7  DEFLECTION  AND  SLOPE  FOR  SEGMENT  2 ',//> 

WRITE! 6. 25) 

25  FORMAT (  1  OX  .  ' ARC  LENGTH'  . 5X ,  ' DEFLECTION  7  . 7X ,  7  SLOPE  7  . 1  OX .  7  STRESS  7 
-  . 1 2X . 7 MOMENT  7  , 1 IX .  7  SHEAR '  .//) 

EXTERNAL  FCN 1 

R= WIDTH- THICK/ ( A  LOG! RA DO /RADI ) ) 

N  W  =  2 

THETA= ARCOS ( (5. *TH I CK- BOND ) / ( 6 . * THICK ) ) 

N*  2 

YBAR=RADI +THICK/2 . - R 
XsO.O 

Y  !  t  )  =  -  T  C  1  ) 

Y  !  2  )  =  -  T  (  2  ) 

TOL=. 000001 

INOsI 

u=o 

DO  510  Z-0.0,R*THETA+.0Q1 , R* THETA/STEP 
XEND=F10AT(Z ) 

CALL  DVE  RK ( N . FCN1 . X . Y , XEND . TOl „ IND.C.NW, W. I.£R ) 

IF! IND. LT.O.OR. IER.GT.O)  GO  TO  5l  1 
dW=d+10 

IF! I NT (JW/ 10). NE . UW/ 1 0.0)  GO  TO  299 
WR I TE( 6 , 30 )  X.Y:i).Y(2).  TORC,  SHEAR 
30  FORMAT! I0X,3(E11.4.3X) .20X .Ell. 4. 5X. Ell. 4,/) 

YGLOB=- . 5*TH1CK 

DO  660  H=-(RAOO-R) .R-RADI .THICK/10. 

----STRESSES  FROM  THE  MOMENT  DISTRIBUTION  SUPERIMPOSED 

- ON  TENSILE  STRESSES  AS  A  FUNCTION  OF  BEAM  THICKNESS 

AREA=WIOTH*THICK 

BSIGX=TORC*H/ ( ( R -H )* WI DTH* THICK* YBAR ) +P LOAD- COS { X/R ) /AREA 


00024900 
00025000 
0002=100 
00025200 
00025300 
00025400 
0  3025500 
03025600 
00025700 
03025300 
03025900 
C  3026000 
03026100 
03026200 
03025220 
0  3026250 
03026251 
03026252 
03026253 
0 3026300 
0  3026-400 
0302650C 
0  30266C0 
C 3026650 
C  3026650 
03026700 
03026800 
03026900 
C 0027000 
03027100 
COO 2 7200 
03027300 
02027400 
00027500 
C 3027600 
C  3027700 
COO 2  7900 
COO 2 8000 
C  3028 1 00 
C 3028200 
C  3028300 
C  3028400 
03026500 
C  3023600 
C 3028700 
C 0028900 
C 0029000 
C  0029 1 00 
C  30  29200 
C 0029200 
C  0029400 
COO 2 9500 
C  0029600 
C  0029700 
C00298C0 
C0029900 
00029540 
C  0029950 
C  0029955 
C0029957 
00029959 
00029960 


o  o  o  o  oo  oooo 
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WR ITE { 6 , 600)  YGL06.  BSIGX 
600  FORMAT (52X.F4.2.1X.E11 .  4 ) 

YGL0B=YGL0 3+ THICK/10 , 

660  CONTINUE 
299  U  =  U+1 

510  CONTINUE 

511  CONTINUE 

- - CALCULATION  OF  DEFLECTION  AND  SLOPE  FOR  SEGMENT  3 

AS  A  FUNCTION  C-  ARC  LENGTH 

WRITE (6. 27) 

27  FORMAT (/////, '  DEFLECTION  AND  SLOPE  FOR  SEGMENT  3*.// I 
WRITE (6. 26) 

26  FORMAT (  1  OX. '  ARC  LENGTH' . 5X .  ' DEFLECTION '  ,7X.  'SLOPE'  . 1 CX , ' STRESS' 
- . 12X. 'MOMENT' .1 IX. 'SHEAR' .//) 

EXTERNAL  FCN2 
X  =  0.0 

Y(l)s-YM) 

Y( 2)«-Y( 2) 

NW  =  2 
N  =  2 
IND«1 

DO  515  M*1 .24 
C(M)aO.O 
515  CONTINUE 
J  =  0 

DO  520  2=0.0. RXTHETA+. 001  . R* THETA /ST E P 

XEND=  FLO AT ( 2 ) 

CALL  DVERK(N. FCN2.X . Y.XEND.TGL. I ND , C , NW . W . I ER ) 

I  F  (  INO.  LT..0.0R.  IER.GT.O)  GO  TO  525 
UW  =  U  +  1  0 

IF( INT( JW/10) .NE.JW/10.0)  GO  TO  349 
WR I TE ( 6 , 60  )  X.  Y(T).  Y(2).  TCRCS.  SHEAR 
60  FORMAT (  1  OX , 3 ( E 1  !  . 4 . 3X  )  . 20X . E  1  1  . 4 . 5X . E 1 1  ,  4 , / ) 

YGLOB=- . 5*TH I CK 

00  770  H=  -  (RADO-R ). R - RADI  . THICK/10 . 

. STRESSES  FROM  THE  MOMENT  DISTRIBUTION  SUPERIMPOSED 

. ON  TENSILE  STRESSES  AS  A  FUNCTION  OF  BEAM  THICKNESS 

CS  IGX=  -  TORCS”  H/  {  (R-H)* WIDTH* THICKLY BAR)  +PLOAD  ’“COS  (THETA-X/R)/ 

-  {  WIDTH ' THICK )  a 

WR I TE ( 6 . 700 )  YGLOB .  CSIGX 
700  FORMAT (52X.F4.2, IX. El  1  .4) 

YGLOB= YGLOB+THI CK/1 0 . 

770  CONTINUE 
349  U  =  U+ 1 
520  CONTINUE 
525  CONTINUE 

. CALCULATION  OF  DEFLECTION  AND  SLOPE  AS  A  FUNCTION 

OF  X  FOR  SEGMENT  4 

WRITE(S, 28) 

28  FORMAT (/////. '  DEFLECTION  AND  SLOPE  FOR  SEGMENT  4',//) 

WR I TE ( 6 . 29) 

29  FORMAT! 1  OX, 'X  DISTANCE 5X DEFLECT  ION 7X SLOPE ',// ) 

EXTERNAL  FCN3 

X  =  0.0 
NW  =  2 
N  =  2 
I  ND=  1 

DO  530  M*1  .24 


0 302996 1 
03029962 
03029966 
03029963 
03030000 
03030100 
03030200 
00030300 
03030400 
0  3030500 
00030600 
03030700 
0  303  0800 
00030900 
03031000 
03031 100 
03031200 
0  303  1300 
00031400 
03031500 
00031700 
03031800 
030319GG 
02032000 
03032100 
00032200 
00032300 
C  303  2400 
C 3032500 
C  303  2600 
03032700 
03032800 
C  0032900 
C  3033000 
00033100 
C  30331 40 
0303315 0 
C  3033 1 55 
C  3  0  3  3 1  5  7 
C  3033 1 60 
C 30331 61 
C 3033 1 62 
C  3  0  3  3  1  6  4 
C  303  3 1 66 
C  3033 1 65 
C3033200 
C 3033300 
C 3033*00 
C 2 03 3600 
C  3033700 
C 003 3800 
C 003 3900 
C 3034000 
00034100 
COO 3 4 200 
C  0034300 
00034400 
COO 3 4500 
C0034600 
C  0034700 
C 3034800 
C  0034900 
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C  (  M )  =  0 . 0 
530  CONTINUE 

vJ  =  0 

DO  525  2=0.0. C0NTA+ . 005 , C0NTA/STEP 
XEND=FL0AT(Z) 

CALL  DVERK(N. FCN3.X . Y .XEND . T0L . IND . C , NW . W . I ER ) 

IFf IND. IT.O.OR. IER.GT.O)  GO  TO  540 
UW=U+10 

I F { INT{ UW/10) .NE. JW/10.0)  GO  TO  399 
WRITE (6. 62)  X.  Y( 1 ) .  Y(2) 

62  FORMAT  (  1  OX  .  3(  E'l  1  . 4 . 3X  )  .  /  ) 

399  u  =  u+: 

535  CONTINUE 
540  CONTINUE 

. CALCULATION  OF  DEFLECTION  AND  SLOPE  AS  A  FUNCTION  OF 

X  FOR  SEGMENT  5 

WRITE (6. 31) 

31  FORMAT {/////, 1  DEFLECTION  AND  SLOPE  FOR  SEGMENT  S',//) 

WRITE (6. 32) 

32  FORMAT {  1  OX . ' X  DISTANCE ' . 5X . ' DEFLECTION '  ,  7X  ,  'SLOPE'  . 1  OX ,  'STRESS' 
-  ,  12X. 'MOMENT' .//) 

EXTERNAL  FCN4 
X  =  0.0 
r;w  =  2 
N  =  2 
I  ND=  1 

DO  545  Mai  .24 
C/M)s0.0 
545  CONTINUE 
J  =  0 

DO  550  Z  =0 . 0  .  SEGB-K 005 .SEGB/STEP 
XEND=FL0AT(2) 

CALL  DVE  RK  ( N . FCN4 . X . Y .XEND . TOL . I  NO . C . NW . W , I ER ) 

IF( IND. LT.O.CR. IER.GT.O)  GO  TO  555 
u  W  =  J  + 1  Q 

I F ( X  .  EC .  .15)  GO  TO  244 
I  F  (  INT(.JW/10)  .NE.  JW/10.0)  GO  TO  245 

244  WR  I  TE  (  6 . 65  )  X.  Yd).  Y(2).  PLGAD'Yd) 

65  FORMAT (  1  OX . 3 ( E 1 1  . 4 , 3X  )  .20X.E11 .4,/) 

DO  880  H  =  -THlCK/2.  . THICK/2  .  .THICK/10. 

. STRESSES  FROM  THE  MOMENT  DISTRIBUTION  SUPERIMPOSED 

. CN  TENSILE  STRESSES  AS  A  FUNCTION  OF  BEAM  THICKNESS 

ESIGXa-PLOAD*YM  )'H/(  1  ./I  2 . "WIDTH* THICK* *3 )  +  P LOAD/ { WIDTH1  THICK  ) 
WR I TE ( 6 . 800 )  H.  CSIGX 
800  FORMAT  (  52X.F4.2,  IX.  Ed  .4) 

880  CONTINUE 

245  J  =  J+ 1 
550  CONTINUE 
555  CONTINUE 

STOP 
.  END 

-4^*.*  *********************  *.*;ir*j»**4#*3jr********#*3'c*»r***ir***tt 

-  ************  ***********. +  *****^****Jjii**)f**r*****  +  ********** 

****  SUBROUTINES  *** 

-**************#*******. -+  +  *^^************3r*^*****J((**  *  **■**,- 

-  ***********  ****#**  +  +  **y#i|f***»^4«***4t^******^#,(ty;jtT^^**!j.jr#^^ 


NUMERICAL  INTEGRATION  OF  SEG2 
SUBROUTINE  FCN1 ( N . X , Y . YPR I  ME ) 


C  9035000 
0 3035100 
03035200 
03035300 
03035400 
03035500 
03035600 
00035700 
03035800 
03035900 
00036000. 
00036100 
03036200 
03036300 
C 0036400 
00036500 
00035600 
00036700 
03026800 
00036900 
C0037CC0 
00037100 
00037150 
00037200 
00U373C0 
C  0037400 
03037500 
00037600 
C  0037700 
0  3037800 
00037900 
00038000 
00038100 
00038200 
C  30383C0 
C0038400 
00038500 
C303S550 
C  0038600 
C 0028700 
C  0038800 
00038840 
C  3036850 
C 0038355 
0 0038060 
00038680 
0 0038890 
C  0038895 
C003S900 
C  0039000 
00039100 
C  3039300 
00039400 
00029410 
00039420 
C0039430 
C  003944Q 
C  0033450 
C  0039460 
COO 3 94 70 
C00394SO 
C  0039500 
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COMMON  R.  P LOAD ,  ESMC  ,  WIDTH.  THICK.  YBAR.  THETA.  EACH.  PI 
BOND.  J.  TORC.  TCRCS ,  TRAC.  TRACE.  SHEAR 
DIMENSION  Y(2).  YPRIME(2) 

YPR IME( t  )=Y{2  } 

-ECCENTRICITY  DUE  ONLY  TO  EXTENSlCNAL  EFFECTS 

YECCsPLOA'D*X*COS(X/R  )  *SIN( X/R )/( WIDTH*  THICK'  E ST/IC  ) 

TORC^PLC AD* ( YECCXTRAC+YBAR+R* ( 1  . -COS( X/R ) ) - Y { 1  )  *COS { X/R ) } 
SHEAR  =  P  LOAD*  {  PLOAD/ ( WIDTH* THICK* ESMC) *  <X/R*COS( X/R  >* *2 
-+S I N ( X/R  )  * ( X / R  * (  -  SI N ( X/R } )+COS(X/R) )  J+SIN( X/R )+Y (  1  ) /R 
-*SIN(X/R)) 

560  YP  R  I  ME  (  2  )«(  R+Y{  t  )  )  *T0RC/ (  R**2  .  +ESMC*WI  DTH'THl  CK*YBAR  ) 

RETURN 

END 

- - NUMERICAL  INTEGRATION  OF  SEG3 

SUBROUTINE  FCN2 (N . X . Y . YPR I  ME ) 

COMMON  R.  PLOAD,  ESMC.  WIDTH,  THICK.  Y BAR ,  THETA,  EACH.  PI 
BOND.  J ,  TORC.  TORCS,  TRAC,  TRACE,  SHEAR 
DIMENSION  Y  (  2.)  ,  YPPIME(2) 

.  YPRIMEf 1  )  =  Y  (  2  ) 

• - ECCENTRICITY  DUE  TO  GEOMETRY 

AEGEO  =  YB AR+R* ( 1  . -COS ( THETA J  )  - 2 . * Y8AR* COS ( THE TA  ) 

. ECCENTRICITY  DUE  TO  GEOMETRY 

3EGEO=  R * ( S I N ( X/R+PI / 2 . -THETA ) -S IN( PI/2 .  -THETA)  ) 

- ECCENTRICITY  DUE  ONLY  TO  EXTENS ION AL  EFFECTS 

EEXT=PLOAD*(SlN( PI/2 . -THETA+X/R ) ) *SIN( THETA- X/R ) *X/( 

-WIDTH '‘THICK- ESMC) 

. ECCENTRICITY  DUE  TO  DEFLECTION 

EDE  FL  =  Y (  1  )*COS*  THETA  -  X /R ) 

TORCS  =  PLOAD* f  AEGEC+BEGLO+E DEFL+EEXT *  TRACE ) 

SHEARaP  LOAD18  (  COS  (X/R  +  P 1/2.  -  THETA)  -  Y  (  1  J  /R  *  S I N  !  P I  /2  . 

-  - THETA  +  X/R  )  +  P LOAD* ( ( SIN<  PI/2 . -THETA+X/R ))  d  SIN{ THETA  - 
- X/R )/(W I DTH*THICK  ‘ ESMC) -X/( WIDTH' THICK*R*ESMC ) ‘COS 

- ( THETA -X/R  )  )+X/( WIDTH* THICK* ESMC' R) * S I N ( THET A  - X/R ) 

-'COS (PI/ 2. -THETA+X/R) ) ) 

Y PRIME (2  )  =  (R+Y{ 1 ) i *  TCRCS/ ( R *  *2 . *  ESMC  * W I DTH* TH ICK*  YBAR ) 

RETURN 

END 

. NUMERICAL  INTEGRATION  OF  SEG4 

SUBROUTINE  FCN3 ( N , X . Y . YPR I  ME ) 

COMMON  R,  PLOAD.  ESMC.  WIDTH.  THICK.  YBAR.  THETA,  EADH .  PI 
- ,  BOND 

DIMENSION  Y(2 ) .  YPRIME(2).  El ( 3 ) 

Eld  ) =ESMC* ( 1 ./I  2. "THICK* *3. *WIDTH+WIDTH*THlCK*{  . 5* TH I CK  +  BOND 

-  /  2  .  )  *  ■*  2  .  ) 

El  (  2)  =  EADH*t ./I  2.* THICK* *3* WIDTH 

El (3)  =  ESMC*( 1  . / 1 2 . *W  I DTH*  THICK*  *3+WIDTHMTHICK* ( . 5* TH I CK+ BOND/ 
-2.  )  *  *  2  .  ) 

YPR IME ( 1  )  aY { 2 ) 


00039600 
00039700 
00039800 
00039900 
00039950 
00039955 
00039960 
00040000 
•00040100 
00040150 
0004Q160 
00040170 
00040500 
o:o4G5oo 
00040700 
00040750 
03040755 
00040760 
03040800 
03040900 
0304 1000 
0304 1 100 
0  3-04  1200 
0304 1 250 
03041260 
0304 1 270 
03041300 
03041350 
0304 1 360 
00041370 
00041400 
00041450 
0004 1 460 
00041470 
00041500 
0004  1600 
C  304 1 650 
C  304 1 660 
0304 1670 
05041700 
C3041 COO 
03041850 
0304 i 850 
05041870 
C 304  1880 
03041830 
05042000 
COO 4 2 1 00 
C  0042200 
C 0042250 
C  0042260 
CD042270 
00042300 
C 0042400 
C  004  2500 
C  004  2600 
C 004 2700 
C  0042800 
C  0042900 
C00430G0 
00043100 
CC043200 


o  o  o 
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DEN  =  EI  (  1  )+EI(2')+EI(3) 

Y  PR  If/IE  (  2  )  UPLOAD*  (  .  5*  THICK+BOND/2  .  +Y  (  1  )  )  /DEN 

RETURN 

END 

-NUMERICAL  INTEGRATION  OF  $EG5 
SUBROUTINE  FCN4{ N . X . Y . YPR I  ME ) 

COMMON  R.  PLOAD,  ESMC .  WIDTH.  THICK.  YBAR.  THETA,  EADH.  PI 
- .  BOND 

DIMENSION  Y ( 2  )  .  YPR I  ME { 2 ) 

YPR IME {  1  ) =Y ( 2  } 

Y PRIME! 2) UPLOAD* (Y{ 1 ) )/(ES MC* WIDTH* THICK* *3. /12. ) 

RETURN 

END 


0 3043300 
Q30434C0 
0  3043500 
0  3043600 
0 D043650 
03043560 
03043670 
03043700 
03043300 
03043300 
03044000 
C 3044100 
03044200 
03044300 
03044400 
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b .  CONVERT 

The  program  CONVERT  essentially  performs  the 
tedious  calculations  involved  in  computing  the  boundary 
conditions  for  the  finite-element  model.  Stresses  dictated 
by  the  beam  bending  model  are  converted  to  equivalent  point 
forces  which  are  then  applied  to  the  finely  meshed  ends  of 
the  finite-element  structure.  In  converting  the  stress 
distribution  from  deformed  to  undeformed  geometry  the 
program  insures  that  the  model  be  maintained  in  equilibrium 
through  the  introduction  of  a  correcting  moment. 

The  important  parameters  utilized  in  the  routine 
are  defined  in  the  nomenclature  section  of  the  program. 
Frequent  comment  cards  are  intended  to  assist  the  user 
in  the  utilization  of  the  program. 
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c- 

C-  STRESS  TRANSFORMATION  PROGRAM 

C  - 

c-  DEVELOPED  BY:  RICHARD  C.  GlVLER 

UNIVERSITY  OF  DELAWARE 
OCT  78  -  JAN  79 


T  FREE 

DIMENSION  20M{5) .  RES{5).  PART ( 4 ) ,  Z0MX(5J,  RESB(5) 

RESX ( 5  )  .  RESY ( 5 ) .  ZOMA ( 5 ) ,  Z0MB(5).  BM0M(5),  ZOMXB(S) 
RESXB(S),  RESYBC5) 

-PARAMETERS  AND  NOMENCLATURE 

--MATERIAL  THICKNESS 
TH I CK= . 1 

C- . MATERIAL  WIDTH 

WI DTH= 1 . 

C . ROTATION  OF  LEFT  HAND  FACE  FROM  UNCEFORMED  GEOMETRY  (RAD) 

DUDSA-- . 01642 

C . ROTATION  OF  RIGHT  HAND  FACE  FROM  UNDEFORMED  GEOMETRY  (RAD) 

DUDSB=. 008420 

C . TOTAL  MOMENT  ON  LEFT  HANO  FACE  (IN.  LBS.) 

TOTMOA= 10.64 

C . SHEAR  ON  LEFT  HAND  FACE 

SH  EARA= 132.6 

C - TOTAL  MOMENT  ON  RIGHT  HAND  FACE  (IN.  LBS.) 

TOTMOB= -  1  .368 
AREA«WIDTH*THiCK 

C . OUTSIDE  RADIUS  IN  INCHES 

RA  DO  =  3 . 5  *TH I CK 

C . INSIDE  RADIUS  IN  INCHES 

RADI =2 . 5  * TH I CK 

C- . RADIUS  OF  CURVATURE  OF  CURVED  MEMBERS 

R* WIDTH* THICK/ ( A  LOG ( RADO/RADI ) ) 

C - --DIFFERENCE  BETWEEN  NEUTRAL  AXIS  AND  CENTROI DAL  AXIS  (IN.) 

YBAR=RAO I+THICK/2 . -R 

C . ADHESIVE  BOND  THICKNESS  IN  INCHES 

BOND= . 03 

C -  -  -  -  -  -  ANGLE  SUBTENDED  BY  CURVED  MEMBERS  IN  RADIANS 
THETA=ARCOS ( ( 5 , * THICK -BOND ) / ( 6 . *THICK ) ) 

Ct--- --TENSILE  LOAD 
PL0AD=200. 

C- 


00000004 
0C000005 
00000006 
00000007 
00000008 
00000010 
0000001 1 
00000012 
COOOOO  t  3 
00000014 
00000015 
00000016 
00000017 
00000018 
00000019 
C0000020 
00000030 
C000003 1 
00000032 
00000033 
00000034 
03000035 
00000036 
0  0000037 
03000100 
00000200 
00000300 
00000400 
00000420 
00000430 
00000440 
C  000047  0 
00000500 
C0000550 
00000600 
00000650 
C 0000700 
C0000750 
00000800 
00000850 
00000900 
00000925 
00000950 
00000975 
00001000 
00001 100 
00001150 
.03001200 
C3001250 
00001300 
C3001 350 
03001400 
00001450 
C0001500 
00001550 
C0001600 
00001650 
00001700 
00001750 
C0001800 
C0001850 


ooooooono  onooo  non  %  oooooonoooo  o  ooo 
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-NOTE:  DEFLECTIONS  ARE  MEASURED  NORMAL  TO  THE  NEUTRAL  AXIS 

-DEFLECTION  OF  LEFT  HAND  FACE  (IN.) 

DEFLA=-. 01547 

-DEFLECTION  OF  RIGHT  HAND  FACE  (IN.) 

DEFLB*- .007872  ' 


-  ************  ********************************************** 


RESOLVING  STRESS  DISTRIBUTION  ON  LEFT  HAND  FACE 


.>**********>********************************************** 


calculation  of  resultant  point  forces  from  the 
STRESS  DISTRIBUTION 


H  =  R -  RAD  I 
DO  100  N  =  1  ,5 

PART {  1  ) = TOTMOA/ ( YBAR * AREA l * ( R-H- R*AL0G( R-H)  ) 

-  +  P LOAD/ A  REA* COS ( THE T  A+DUD5 A ) *H 
HsH- . 02 

PART(2)sT0TM0A/(YBAR*AREA)*(R-H-R*AL06(  R-H)  ) 

-  +  P  LOAD/ A  REA* COS ( THET A+DUDS A ) *H 
RES l  N  )  =  PART ( 1  ) -  PART ( 2 ) 

100  CONTINUE 


CALCULATION  OF  ACTUAL  MOMENTS  FROM  THE  STRESS  DISTRIBUTION 


h=r-radi 

DO  200  N  =  1  .5 

PART<3)=T0TM0A/(Y8AR*AREA)*(  -1  .  )M  .  5M  R -H  )  *  *2  -  2  .  *R*  <  R-H  )+R-*"2 
-  *  A  LOG  (  R-  H)  )  +  .  5* P LOAD /A RE A' COS (  THETA+DUDSA)*H**2 
H  =  H -  .02 

PART (4 )r TOTMOA/( YBAR * AREA )*(-1.)+(.5*(R-H)**2-2.*R* ( R-H )+R**2 
-*  A  LOG( R  -  H ) )+ . 5*  P LOAD/ AREA*  COS ( THETA+DUDSA )  *2 

ZO.M  {  N  )  =  P  ART  (  3  }  -  PART  f  4  ) 

200  CONTINUE 


CALCULATION  OF  CORRECTION  MOMENT  DUE  TO  THE 
REPRESENTATION  OF  THE  STRESS  DISTRIBUTION  BY  POINT 
FORCES 


H=: R  -RADI  -  .01 
DO  300  N  =  1  .5 
Z0MX(N)=ZQM<N)-RES(N)*H 
H=H - .02 
300  CONTINUE 


CALCULATION  OF  MOMENT  DUE  TO  TRANSLATION  OF  THE  STRESS 
DISTRIBUTION  THROUGH  SPACE  FROM  THE  DEFORMED  GEOMETRY 
TO  THE  UNDEFORMED  GEOMETRY 


CALCULATION  OF  MOMX  FOR  INPUT  INTO  THE  FINITE  ELEMENT 
MODEL 


H=  R -RAD I  - .01 
00  400  N  *  1 , 5 

ZOMA{N)=-RES(N)'*(  (H  +  OEFLA)  *C0S(DUDSA)  -  H) 
ZOMX ( N ) * ZOMX ( N ) +ZOMA ( N ) 

HaH- .02 


00001852 
00001854 
00001870 
00001900 
00001970 
00002000 
00002010 
00002020 
00002100 
00002200 
00002300 
00002320 
00002400 
00002500 
00002600 
00002700 
00002800 
00002900 
00003000 
00003100 
00003200 
00003300 
00003400 
C  0003500 
00003600 
00003800 
00003900 
00004000 
00004100 
00004200 
00004300 
00004400 
00004600 
00004700 
C  0004800 
00004900 
00005000 
00005200 
00005300 
00005400 
00005500 
00005600 
00005700 
C0005800 
00005900 
00006000 
00006200 
00006300 
C0006400 
00006500 
C0006600 
00006700 
CC006800 
00006900 
C0007000 
00007100 
00007200 
C0007300 
00007400 
00007500 
CC0077Q0 
00007900 
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400  CONTINUE 


CALCULATION  OF  FY  ANO  FZ  FOR  INPUT  INTO  THE  FINITE 
ELEMENT  MODEL 


WRITE(6. 25) 

25  FORMAT!///////, 6X, 'BOUNDARY  CONDITIONS  FOR  LEFT  HAND  SEGMENT 
WRITE (6. 30) 

30  FORMAT!/  ,4X, 'NODE'  ,6X,  '  FY '  , 1 3X . 'FZ' .12X. 'MOMX'  , // ) 

DO  500  N*1 .5 

RESX ( N )  =  -  RES ( N ) *C0S ( THET A+OUDSA ) -  SHEAR A/5 . *SI N ( THETA ) 

RESY ( N )  =  -  RES! N ) *SIN ( THETA  +  DUDSA ) +SHEARA/5. *COS ( THETA  ) 
WRITEIS.2)  N.  RESX(N).  RESY(N),  - ZOMX ( N ) 

2  FORMAT ( 5X . F2 . 0 , 3X . 3! El 0. 4 . 5X ) , / ) 

500  CONTINUE 


•  *■*****+*****  ******* |r  *  *  *  *  *  ************************.****■**  +  **;* 


RESOLVING  STRESS  DISTRIBUTION  ON  RIGHT  HAND  FACE 


.*«*********************»*********************  ************ 


CALCULATION  OF  RESULTANT  POINT  FORCES  FROM  THE  STRESS 
DISTRIBUTION 


H=  TH I CK/ 2 . 

DO  1 000  N= 1 . 5 

PART  (  1  )  =  *  .5*T0TM0B*H**2/(  1  .  / 1  2  .  *  W I DTH*  TH  I  CK  **  3  )  +  P  LOAD/ AR  E  A  *H 
H=H - .02 

PART ( 2 )  =  - . 5*T0TM08*H*  *2/ (  1  . /I  2.  * W I DTK* TH I CK *  *  3 ) +P LOAD/ AREA*H 
RESB ( N )  =  PART (  1  ) - PART ( 2 ) 

1000  CONTINUE 


CALCULATION  OF  ACTUAL  MOMENTS  FROM  THE  STRESS 
DISTRIBUTION 


H» THICK/2. 

DO  1100  N= 1 , 5  o 

PART ( 3 )  =  1  . /3 . *(T0TM0B)*H**3/(1 . / 1 2 . *WI DTH*TH ICK*  *  3 ) 
-  - . 5*PL0AD/AREA*H**2 
HsH- .02 

PART(4)  =  1 ./3. M TOTMOB ) *H*  *  3/ { 1 . /I  2 . *WI DTH*THICK**3 ) 
-- . 5*PL0AD/AREA*H**2 
BMOM (  N  )  =  PART ( 3 ) -  PART ( 4 ) 

1  100  CONTINUE 


CALCULATION  OF  CORRECTION  MOMENT  DUE  TO  THE 
REPRESENTATION  OF  THE  STRESS  DISTRIBUTION  BY  POINT 
FORCES 


H«THICK/2. - .01 
DO  1200  N= 1 .5 

ZOMXB ( N ) =  BMOM (N)+RES8(N)*H 
H  =  H -  .02 
1200  CONTINUE 


CALCULATION  OF  THE  MOMENT  DUE  TO  TRANSLATING  THE  STRESS 
DISTRIBUTION  FROM  DEFORMED  TO  UNDEFORMED  GEOMETRY 


00008000 
00008100 
00003200 
00008300 
00008400 
00008500 
/)  00008600 
00008700 
00008800 
00008900 
00009000 
00009100 
00009200 
00009300 
00009400 
00009500 
00009600 
00009700 
00009800 
00009900 
00010000 
00010100 
00010200 
00010300 
00010400 
00010500 
00010600 
00010700 
00010800 
00010900 
0001 1000 
0001 1 100 
0001 1200 
0001 1300 
0001 1 500 
0001 1600 
0001 1700 
0001 1800 
0001 1900 
00012000 
00012100 
CO0 1 2200 
00012300 
C001 2350 
00012400 
00012500 
C001 2600 
C001 2800 
00012900 
CO0 1 3000 
C001 3100 
00013200 
C001 3300 
C001 3400 
00013500 
CO0 1 3600 
00013750 
00013800 
0001 3900 
C0014000 
COO 14100 
CO0 1 4200 
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CALCULATION  OF  MOMENT  MX  FOR  INPUT  INTO  THE  FINITE 
ELEMENT  MODEL 


H  =  THICK/ 2 . - .01 
DO  1300  N=  1  » 5 

ZOMB(N)=RESB(N)*( ( H- DEFLB ) *C0S( DUDSB ) - H ) 
ZOMXB (  N  )■  »Z0MXB(  N  )+Z0M8( N  ) 

H=H- .02 
1300  CONTINUE 


CALCULATION  OF  FY  AND  FZ  FOR  INPUT  INTO  THE  FINITE 
ELEMENT  MODEL 


WRITE (6. 40) 

40  FORMAT (//////, 1  OX . ' BOND ARY  CONDITIONS  FOR  SEG8',/) 
WR I TE ( 6 , 45  ) 

^45  FORMAT  (  /  .4X,  '  NODE'  ,6X,  '  FY#  .  1  3X  .  '  FZM2X,  'MOMX'  ,  //  ) 
DO  1400  N=1 ,5 

RESXB ( N ) =RESB( N ) *C0S { DUDSB ) 

RESYB(N) =RESB(N)*SIN(DUDS8 ) 

WR  I  TE ( 6 , 20 )  N.  RESXB ( N ) ,  RESYB(N).  ZOMXB(N) 

20  FORMATS  5X,F2.0,3X,3(E10.4.5X),/) 

1400  CONTINUE 
END 


00014300 
00014400 
000145G0 
00014600 
00014700 
00014800 
00014900 
00015100 
00015200 
00015300 
00015400 
00015500 
00015600 
00015700 
00015800 
00015900 
00016000 
00016100 
00016200 
COO 1 6300 
00016400 
00016500 
00016600 
00016700 
00016800 
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DVERK 


N 

FCN 


X 


Y 


XEND 


TOL 


(N . FCN , X . Y .XEND, TOL.  I  NO , C . NW , W . I ER } 

LIBRARY  3 . . . 

-  SOLUTION  OF  A  SYSTEM  OF  FIRST  ORDER  ORDINARY 
DIFFERENTIAL  EQUATIONS  OF  THE  FORM 
DY/DX  =  F { X . Y )  WITH  INITIAL  CONDITIONS. 

A  runge-kutta  METHOD  BASED  ON  verners  fifth 
AND  SIXTH  ORDER  PAIR  OF  FORMULAS  IS  USED. 

*  CALL  DVERK(N. FCN. X.Y, XEND. TOL, IND. C.NW.W.IER) 

*  NUMBER  OF  EQUATIONS.  (INPUT) 

*  NAME  OF  SUBROUTINE  FOR  EVALUATING  FUNCTIONS. 

( INPUT  ) 

THE  SUBROUTINE  ITSELF  MUST  ALSO  SE  PROVIDED 
BY  THE  USER  AND  IT  SHOULD  BE  OF  THE 
FOLLOWING  FORM 

SUBROUTINE  FCN(N.X. Y.YPRIME) 

DIMENSION  YIN) .YPRIME(N) 


FCN  SHOULD  EVALUATE  YPRIME(I) . YPRIWE(N) 

GIVEN  N.X.  AND  Y(t) . Y ( N  )  .  YPRIME(I) 

IS  THE  FIRST  DERIVATIVE  OF  Y ( I )  WITH 
RESPECT  TO  X. 

FCN  MUST  APPEAR  IN  AN  EXTERNAL  STATEMENT  IN 

THE  CALLING  PROGRAM  AND  N.X.Y(I) . Y(N) 

MUST  NOT  BE  ALTERED  BY  FCN. 

-  INDEPENDENT  VARIABLE.  (INPUT  AND  OUTPUT) 

ON  INPUT.  X  SUPPLIES  THE  INITIAL  VALUE. 

ON  OUTPUT.  X  IS  REPLACED  WITH  XEND  UNLESS 
ERROR  CONDITIONS  ARISE.  SEE  THE  DES¬ 
CRIPTION  OF  PARAMETER  IND. 

-  DEPENDENT  VARIABLES.  VECTOR  OF  LENGTH  N. 

( INPUT  AND  OUTPUT  ) 

ON  INPUT,  Y  (  1  ) . Y  (  N  )  SUPPLY  'INI  T  I AL 

VALUES. 

ON  OUTPUT.  Y I  1  ) . Y ( N  )  ARE  REPLACED  WITH 

AN  APPROXIMATE  SOLUTION  AT  XEND  UNLESS 
ERROR  CONDITIONS  ARISE.  SEE  THE  DES¬ 
CRIPTION  OF  PARAMETER  IND. 

-  VALUE  OF  X  AT  WHICH  SOLUTION  IS  DESIRED. 

( INPUT) 

XEND  MAY  BE  LESS  THAN  THE  INITIAL  VALUE  OF 
X  . 

-  TOLERANCE  FQP  ERROR  CONTROL.  (INPUT) 

THE  SUBROUTINE  ATTEMPTS  TO  CONTROL  A  NORM 
OF  THE  LOCAL  ERROR  IN  SUCH  A  WAY  THAT  THE 
GLOBAL  ERROR  IS  PROPORTIONAL  TO  TOL. 

MAKING  TOL  SMALLER  IMPROVES  ACCURACY  AND 
MORE  THAN  ONE  RUN.  WITH  DIFFERENT  VALUES 
OF  TOL.  CAN  BE  USED  IN  AN  ATTEMPT  TO 
ESTIMATE  THE  GLOBAL  ERROR. 

IN  THE  DEFAULT  CASE  (IND=1).  THE  GLOBAL 
ERROR  IS 

MAX ( ABS  (  E  (  1  ) ),  .  .  .  , ABS ( E ( N ) )  ) 

WHERE  E(K)*(Y(K)-YT(K) )/MAX( 1 .ABS(Y(K) ) ) 

YT ( K )  IS  THE  TRUE  SOLUTION.  AND 
Y ( K )  IS  THE  COMPUTED  SOLUTION  AT  XEND. 

FOR  K»1  .2 . N. 

OTHER  ERROR  CONTROL  OPTIONS  ARE  AVAILABLE. 


D/EKOOIO 
D7EK0020 
DV  EK0030 
D7EK0C40 
D/EK0050 
DVEK0060 

d;ekoo7o 
d;ekoo3o 
d;ekooso 
DVEKOIOO 
D^EKOI 10 
DyEKOI 20 
DVEK0130 
DVEK0140 
DVEK0150 
C7EK0160 
D  v'EKOI  70 
0  i/EK01  SO 
D7EK0190 
0VEK0200 
D  */EK02  1  0 
C yEK0220 
C*/EK0230 
Cv/EK0240 
C  n/  EK0250 
C  y/EK0260 
CVEK0270 
D  y/EK0280 
Dn/EKG2S0 
D  y£K0300 
C  s/  EK03  1  0 
DVEK0320 
CVEK0330 
C v EK0340 
C  y/EK0350 
DVEK0360 
Dv/EK0370 
C  EKC380 
CVEK0390 
D  v EK0400 

D</EK041  0 

CVEKC420 
C  i/ EK0430 
C  M EK0440 
CVEK0450 
CVEK0460 
Z VEKC470 
CVEKC4S0 
CVEK0490 
DVEK0500 
0VEKC510 
0  VEK  0520 
DVEK0530 
CVEK0540 
CVEK0550 
CVEK0560 
CVEK0570 
CVEK0580 
CVEK0590 
CVEK0600 
CVEK0610 
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SEE  THE  DESCRIPTION  OF  PARAMETERS  IND  AND 
C  BELOW . 

IND  -  INDICATOR.  (INPUT  AND  OUTPUT) 

ON  INITIAL  ENTRY  IND  MUST  BE  SET  EQUAL  TO 
EITHER  1  OR  2. 

IND  =  1  CAUSES  ALL  DEFAULT  OPTIONS  TO  BE 
USED  AND  ELIMINATES  THE  NEED  TO  SET 
SPECIFIC  VALUES  IN  THE  COMMUNICATIONS 
VECTOR  C- 

IND  =  2  ALLOWS  OPTIONS  TO  BE  SELECTED.  IN 
THIS  CASE.  THE  FIRST  9  COMPONENTS  OF  C 
MUST  BE  INITIALIZED  TO  SELECT  OPTIONS  AS 
DESCRIBED  BELOW. 

THE  SUBROUTINE  WILL  NORMALLY  RETURN  WITH 

IND  =  3.  HAVING  REPLACED  THE  INITIAL  VALUES 
OF  X  AND  Y  WITH.  RESPECTIVELY.  THE  VALUE 
XEND  AND  AN  APPROXIMATION  TO  Y  AT  XEND. 

THE  SUBROUTINE  CAN  BE  CALLED  REPEATEDLY  WITH 
NEW  VALUES  OF  XEND  WITHOUT  CHANGING  ANY 
OF  THE  OTHER  PARAMETERS. 

THREE  ERROR  RETURNS  ARE  ALSO  POSSIBLE.  IN 
WHICH  CASE  X  AND  Y  WILL  BE  THE  MOST 
RECENTLY  ACCEPTED  VALUES. 

IND  =  -3  INDICATES  THAT  THE  SUBROUTINE  WAS 
UNABLE  TO  SATISFY  THE  ERROR  REQUIREMENT. 
THIS  MAY  MEAN  THAT  TOL  IS  TOO  SMALL. 

IND  =  -2  INDICATES  THAT  THE  VALUE  OF  HMIN 
(MINIMUM  STEP-SIZE)  IS  GREATER  THAN  HM AX 
(MAXIMUM  STEP-SIZE).  WHICH  PROBAELY  MEANS 
THAT  THE  REQUESTED  TOL  (WHICH  IS  USED  IN 
THE  CALCULATION  OF  KMIN)  IS  TOO  SMALL. 

IND  =  -1  INDICATES  THAT  THE  ALLOWED  MAXIMUM 
NUMBER  OF  FCN  EVALUATIONS  HAS  BEEN 
EXCEEDED.  THIS  CAN  ONLY  OCCUR  IF  OPTION 
C ( 7 )  .  AS  DESCRIBED  BELOW.  HAS  BEEN  USED. 

C  -  COMMUNICATIONS  VECTOR  OF  LENGTH  24.  (INPUT  IF 

IND. NE. 1 .  AND  OUTPUT ) . 

C  IS  USED  TO  SELECT  OPTIONS  AND  TO  RETAIN 
INFORMATION  BETWEEN  CALLS.  THE  USER  NEED 
NOT  BE  CONCERNED  WITH  THE  FOLLOWING 
DESCRIPTION  OF  THE  ELEMENTS  OF  C  WHEN 
DEFAULT  OPTIONS  ARE  USED  (IND=1). 

HOWEVER.  WHEN  IT  IS  DESIRED  TO  USE  IND-2 
AND  SELECT  OPTIONS,  A  BASIC  UNDERSTANDING 
OF  DVERK  IS  REQUIRED.  THE  FOLLOWING 
PARAGRAPH  DESCRIBES.  BRIEFLY.  THE  BASIC 
TERMS.  FOR  MORE  DETAILS.  SEE  THE 
REFERENCE  - 

DVERK  ADVANCES  THE  INDEPENDENT  VARIABLE 
X  ONE  STEP  AT  A  TIME  UNTIL  XEND  IS 
REACHED.  THE  SOLUTION  IS  COMPUTED  AT 
XTRIAL  =  X+HTRJAL  ALONG  WITH  AN  ERROR 
ESTIMATE  EST.  IF  EST  IS  LESS  THAN  OR 
EQUAL  TO  TOL  (SUCCESSFUL  STEP),  THE  STEP 
IS-  ACCEPTED  And  X  IS  ADVANCED  TO  XTRUL. 
IF  EST  IS  GREATER  THAN  TOL  (FAILURE) 
HTRIAL  IS  ADJUSTED  AND  THE  SOLUTION  IS  . 
RECOMPUTED.  HMAG  =  ABS (HTRIAL)  IS  NEVER 
ALLOWED  TO  EXCEED  HMAX  NOR  IS  IT  ALLOWED 
TO  BECOME  SMALLER  THAN  HMIN.  THE  FIRST 
TRIAL  STEP  IS  HSTART .  DURING  THE 
COMPUTATION.  A  COUNTER  (0(23))  IS 


07 £K 0620 
DVEK0630 
DYEK0640 
D/EKCG50 
D/EK0660 
D/EK0670 
d;eko6so 

D;EK0690 
D/EK0700 
D/EK0710 
D/EK0720 
DVEK0730 
D  t/EK0740 
DVEKG750 
D;EK0760 
D;EK0770 
D  v EK0780 
D  s/EK 0790 
C./EKC800 
D  j  EK08 1 0 
D  v'  EK0820 
Ds/EK.0330 
Dv'EK0840 
Di/EK0S50 
Dv'EKC860 
D  </ EK  0370 
0  7EK0S80 
D  V  EK0330 
DVEK0900 
DVEK0910 
C  v1’  EK  0920 
DVEK0930 
DVEK0940 
Cv'EKCSSO 
CVEKC960 
C^EK0970 
DYEK0980 
D  \/EKC990 
DVEK 1C00 
CVEK’IOIO 
DVEK 1020 
CVEK1030 
CVEK1040 
DVEK 1050 
DVEK 1060 
C  V  E  K  1  0  7  0 
CvEK 1 080 
CVEK1090 
CVEK1 100 
CVEK1 1 10 
DVEK  1  1 20 
CVEK1 130 
CVEK1 140 
DVEK 11 50 
CVEK1 160 
CVEK1 1 70 
CVEK1 ISO 
CVEKI 190 
DVEK 1 200 
DVEK1210 
DVEK 1 220 
DVEK 1 230 


-111- 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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INCREMENTED  EACH  TIME  A  TRIAL  STEP  FAILS  D/EK1240 
TO  PROVIDE  A  SOLUTION  SATISFYING  THE  E R RORD /E K 1 250 


TOLERANCE.  ANOTHER  COUNTER  ( C ( 22 )  )  IS  D/EK1260 
USED  TO  RECORD  THE  NUMBER  OF  SUCCESSFUL  D;EK1270 
STEPS.  AFTER  A  SUCCESSFUL  STEP.  C(23)  IS  D7EK1280 
SET  TO  ZERO.  D/EK1290 

OPTIONS.  IF  THE  SUBROUTINE  IS  ENTERED  WITH  D</EK1300 
I ND=  2  ,  THE  FIRST  9  COMPONENTS  OF  THE  D^EK1310 

COMMUNICATIONS  VECTOR  MUST  BE  INITIALIZED  DVEK1320 
BY  THE  USER.  NORMALLY  THIS  IS  DONE  BY  DVEK1330 

FIRST  SETTING  THEM  ALL  TO  ZERO.  AND  THEN  D,/EK1340 

THOSE  CORRESPONDING  TO  PARTICULAR  OPTIONS  DyEKl350 
ARE  MADE  NON-ZERO.  DVEK 1360 

CM)  -  ERROR  CONTROL  INDICATOR.  D;EK1370 

THE  SUBROUTINE  ATTEMPTS  TO  CONTROL  A  NORM  DVEK1330 

OF  THE  LOCAL  ERROR  IN  SUCH  A  WAY  THAT  THE  DyEKl390 

GLOBAL  ERROR  IS  PROPORTIONAL  TO  TOL.  DVEK1400 

THE  DEFINITION  OF  GLOBAL  ERROR  FOR  THE  DVEK1410 

DEFAULT  CASE  { I N D  = 1 )  IS  GIVEN  IN  THE  D^EK1420 

DESCRIPTION  OF  PARAMETER  TOL.  THE  DEFAULT  DVEK1430 
WEIGHTS  ARE  1  /MAX  (  1  .  A8S  (  Y  (  K  )  )  )  .  WHEN  IND=2  DVEKM40 
IS  USED.  THE  WEIGHTS  ARE  DETERMINED  D VEK 1 450 

ACCORDING  TO  THE  VALUE  OF  CM).  DVEK1460 

IF  C  (  1  )*  1  THE  YiEIGHTS  ARE  1  D7EK1470 

(ABSOLUTE  ERROR  CONTROL)  DvEK 1 480 

IF  C  (  1  )  =2  THE  WEIGHTS  ARE  1/ABS(Y(K))  D^/EK1490 

FOR  K=1  .2 . N  .  DVEK  1  SCO 

(RELATIVE  ERROR  CONTROL)  DVEK1510 

IF  C  (  1  )  =  3  THE  WEIGHTS  ARE  DVEK1520 

1 /MAX ( A BS ( C ( 2  )  )  . A 3 S ( Y ( K  }  )  )  D  VEK  1 530 

.  FOR  K  =  1  ,  2 . N  .  D  EK  1  540 

(RELATIVE  ERROR  CONTROL.  UNLESS  DVEK 1550 

ABS ( Y ( K ) )  IS  LESS  THAN  THE  FLOOR  CWEK1560 
VALUE. A8S(C(2) ) )  DVEK1570 

IF  C ( 1  )  =4  THE  WEIGHTS  ARE  CV EK1580 

1 /MAX ( ABS ( C  f  K+30 )  )  , AES  f  Y ( K  )  )  )  DVEK) 590 

FOR  K  =  1  ,  2 . N.  *  0 VEK 1 500 

(HERE  INDIVIDUAL  FLOOR  VALUES  DVEK1610 

ARE  USED)  C VEK 1 620 

IN  THIS  CASE,  THE  DIMENSION  OF  C  DVEK1630 
MUST  BE  GREATER  THAN  OR  EQUAL  TO  CVEK1640 

N  +  30  AND  C ( 3 1  )  .  C(32) . C1N+30)  CVEK1650 

MUST  be  INITIALIZED  5Y  THE  USER.  CVEK1660 
IF  C(1)»5  THE  WEIGHTS  ARE  1 / ABS ( C ( K  +  30 )  )  C V E K 1 6 7 0 

FOR  K=1  .2 . N.  CVEK16S0 

IN  THIS  CASE.  THE  DIMENSION  OF  C  DVEK1690 
MUST  BE  GREATER  THAN  OR  EQUAL  TO  CVEK1700 

NK*0  AND  C  (  3  1  )  .  C  (  32  > . ClN+30)  CVEK1710 

MUST  BE  INITIALIZED  BY  THE  USER.  CVEK1720 
FOR  ALL  OTHER  VALUES  OF  CM).  INCLUDING  DVEK1730 

C ( 1  )*0  THE  DEFAULT  VALUES  OF  CVEK1740 

THE  WEIGHTS  ARE  TAKEN  TO  BE  C-VEK1750 

t /MAX { 1 . ABS ( Y ( K ) ) )  CVEK1760 

FOR  K=  1  ,  2 . N  .  DVEK  I  770 

0(2)  -  FLOOR  VALUE.  USED  WHEN  THE  INDICATOR . CM )  CVEK1780 

HAS  THE  VALUE  3.  CVEK1790 

C(3)  -  HMIN  SPECIFICATION.  IF  NOT  ZERO.  THE  SUB-  DVEK1800 

ROUTINE  CHOOSES  HMIN  TO  BE  ABS(C(3)).  CVEK1810 

OTHERWISE  IT  USES  THE  DEFAULT  VALUE  CVEK1820 

10-MAX (DWARF , RRE B*MAX ( NORM ( Y ) /TOL . ABS ( X ) )  )  DVEK 1830 
WHERE  DWARF  IS  a  VERY  SMALL  POSITIVE  MACH  I NEC VEK 1 840 
NUMBER  AND  RREB  IS  THE  RELATIVE  ROUNDOFF  CVEK1850 
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CRROR  8OUN0. 

C ( 4 )  -  HST ART  SPECIFICATION.  IF  NOT  ZERO.  THE  SUB¬ 

ROUTINE  WILL  USE  AN  INITIAL  HMAG  EQUAL  TO 
ABS ( C ( 4 ) ) .  EXCEPT  OF  COURSE  FOR  THE  RE¬ 
STRICTIONS  IMPOSED  BY  HMIN  AND  HMAX. 
OTHERWISE  IT  USES  THE  DEFAULT  VALUE 
HMAXMTOL)^  (  1/6)  . 

C ( 5  )  *  SCALE  SPECIFICATION.  THIS  IS  INTENDED  TO  BE 

A  MEASURE  OF  THE  SCALE  OF  THE  PROBLEM. 
LARGER  VALUES  OF  SCALE  TEND  TO  MAKE  THE 
METHOD  MORE  RELIABLE.  FIRST  BY  POSSIBLY  RE¬ 
STRICT  I NG  HMAX  (AS  DESCRIBED  BELOW)  AND 
SECOND.  BY  TIGHTENING  THE  ACCEPTANCE 
REQUIREMENT.  IF  C ( 5 )  IS  ZERO.  A  DEFAULT 
VALUE  OF  1  IS  USED.  FOR  LINEAR  HOMOGENEOUS 
PROBLEMS  WITH  CONSTANT  COEFFICIENTS.  AN 
APPROPRIATE  VALUE  FOR  SCALE  IS  A  NORM  OF 
THE  ASSOCIATED  MATRIX.  FOR  OTHER  PROBLEMS, 
AN  APPROXIMATION  TO  AN  AVERAGE  VALUE  OF'  A 
NORM  OF  THE  JACOBIAN  ALONG  THE  TRAJEC¬ 
TORY  MAY  BE  APPROPRIATE. 

C ( 6 )  '  -  HMAX  SPECIFICATION.  FOUR  CASES  ARE  POSSIBLE, 

IF  C ( 6  )  . NE . 0  AND  C(5).NE.O.  HMAX  IS  TAKEN 
TO  BE  M I N ( ABS ( C ( 6 )  )  . 2/ABS  f  C ( 5  )  )  )  . 

IF  C ( 6  )  . N £ . 0  AND  C ( 5 )  . EO . 0 .  HMAX  IS  TAKEN 
TO  BE  A 8 S  t  C  f  6  )  )  . 

IF  C ( 6  )  . EO . 0  AND  C ( 5 )  . NE . 0 .  HMAX  IS  TAKEN 
TO  BE  2 / A BS ( C ( 5)  )  . 

IF  C ( 6  )  . EQ ■ 0  AND  C{5).EQ.O.  HMAX  IS  GIVEN 
A  DEFAULT  VALUE  OF  2 . 

C ( 7 )  -  MAXIMUM  NUMBER  OF  FUNCTION  EVALUATIONS.  IF 

NOT  ZERO.  AN  ERROR  RETURN  WITH  I NO  =  -1 
WILL  BE  CAUSED  WHEN  THE  NUMBER  OF  FUNCTION 
EVALUATIONS  EXCEEDS  ABS(C(7)). 

C ( 8 )  -  INTERRUPT  NUMBER  1  .  IF  NOT  ZERO.  THE  SUB¬ 

ROUTINE  WILL  INTERRUPT  THE  CALCULATIONS 
AFTER  IT  HAS  CHOSEN  ITS  PRELIMINARY  VALUE 
OF  HMAG.  AND  JUST  BEFORE  CHOOSING  HTR I A  L 
AND  XTRIAL  IN  PREPARATION  FOR  TAKING  A  STEP 
(HTRIAL  MAY  DIFFER  FROM  HMAG  IN  SIGN.  AND 
MAY  REQUIRE  ADJUSTMENT  IF  XEND  IS  NEAR). 

THE  SUBROUTINE  RETURNS  WITH  IND  =  4,  AMD 
WILL  RESUME  CALCULATION  AT  THE  POINT  OF 
INTERRUPTION  IF  RE-ENTERED  WITH  IND  =  4. 

C( 9 )  -  INTERRUPT  NUMBER  2.  IF  NOT  ZERO.  THE  SUB¬ 

ROUTINE  WILL  INTERRUPT  THE  CALCULATIONS 
IMMEDIATELY  AFTER  IT  HAS  DECIDED  WHETHER  OR 
NOT  TO  ACCEPT  THE  RESULT  OF  THE  P^OST  RECENT 
TRIAL  STEP.  WITH  IND  =  5  IF  IT  PLANS  TO 
ACCEPT.  OR  IND  =  6  IF  IT  PLANS  TO  REJECT. 
Y(+)  IS  THE  PREVIOUSLY  ACCEPTED  RESULT. 
WHILE  W  (  *  ,  9  )  IS  THE  NEWLY  COMPUTED  TRIAL 
VALUE.  AND  W ( * . 2 )  IS  THE  UNWEIGHTED  ERROR 
ESTIMATE  VECTOR.  THE  SUBROUTINE  WILL  RESUPvlE 
CALCULATIONS  AT  THE  POINT  OF  INTERRUPTION 
ON  RE-ENTRY  WITH  IND  =  5  OR  6. 

IND  MAY  BE  CHANGED  BY  THE  USER  IN  ORDER  TO 
FORCE  ACCEPTANCE  OF  A  STEP  (BY  CHANGING  IND 
FROM  6  TO  5)  THAT  WOULD  OTHERWISE  BE 
REJECTED.  OR  VICE  VERSA. 

NW  -  THc  FIRST  DIMENSION  OF  W  AS  IT  APPEARS  IN  THE 

CALLING  PROGRAM.  (INPUT) 


D  y  EK  1 860 
D/EK1870 
DVEK1880 
D  y  EK  1 890 
DVEK 1900 
DVEK1910 
DVEK1920 
D7EK1930 
DVEK 1940 
D  EK  1  950 
DVEK1960 
D  y  EK 1 970 
DVEK 1930 
DVEK 1 990 
DVEK2000 
D  y EK  20 1 0 
DVEK 2020 
D  y EK2030 
D;EK2040 
D  y  EK2050 
Dy£K2C60 
DVEK2070 
DVEK2C8Q 
D  yEK2090 
D VEK2100 
0VEK21 10 
DVEK2120 
DVEK2130 
DVEK2140 
0VEK2150 
DVEK2160 
DVEK2170 
DVEK21S0 
CVEK21 90 
CVEK2200 
DVEK221 0 
C  yEK2220 
DVEK2230 
D  v EK  2240 
D  y EK2250 
DVEK2260 
CVEK2270 
CVEK2280 
C V EK2230 
CVEK2300 
CVEK2310 
DVEK2320 
DVEK2330 
DVEK2340 
C V  EK  2350 
CVEK2350 
CVEK2370 
CVEK2380 
Cy£K2390 
CVEK2400 
CVEK2410 
CVEK2420 
CVEK2430 
CVEK2440 
CVEK2450 
C  VEK  2460 
CVEK2470 


nooooooooooooooooooooo 


-113- 


NW  MUST  BE  GREATER  THAN  OR  EQUAL  TO  N. 

W  -  WORKSPACE  MATRIX. 

THE  FIRST  DIMENSION  OF  W  MUST  BE  NW  AND  THE 
SECOND  MUST  BE  GREATER  THAN  OR  EQUAL  TO  9. 
IER  -  ERROR  PARAMETER.  (OUTPUT) 

TERMINAL  ERRORS 


IER  = 

129. 

NW  IS  LESS  THAN  N  OR  TOL  IS  LESS 
THAN  OR  EQUAL  TO  ZERO. 

IER  = 

130. 

IND  IS  NOT  IN  THE  RANGE  1  TO  6 . 

IER  = 

131  . 

XEND  HAS  NOT  BEEN  CHANGED  FROM 
PREVIOUS  CALL  OR  X  IS  NOT  SET  TO 
THE  PREVIOUS  XEND  VALUE. 

IER  = 

132. 

THE  RELATIVE  ERROR  CONTROL 

OPTION  ( C ( 1  )  =  2 )  WAS  SELECTED  AND 
ONE  OF  THE  SOLUTION  COMPONENTS 

IS  ZERO. 

PRECISION 

-  SINGLE 

REOD.  IMS L 

ROUTINES  -  UERTST 

LANGUAGE 

*  FORTRAN 

LATEST  REVISION 

.  INTEGER 
INTEGER 
REAL 
REAL 

1 

REAL 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

OATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


DECEMBER  15.  1976 

BGH 

N, IND.NW.K 
IER 

X , Y { N ) . XEND . TOL ,C{ 1  )  . W ( N W . 9  )  . TEMP 

ZERO. ONE . TWO. THREE . FOUR . FIVE . SEVEN . TEN .HALF . P9 

C4D15.C2D3.C5D6.C1D6.C1D1 5 .C2D96 

RK ( 39) . REPS. RTOL 

ZERO/O.  0/  .ONE/1  .0/  ,  TWO/2 . 0/  .  THRE.E/3 .0/ 
FOUR/4.0/. FI  VE/5.0/. SEVEN/7 ,0/ 

TEN/10.0/ .HALF/0.5/ , P9/0 . 9/ 

C4D 1 5/ . 26666566657/ 

C2D3/. 66666666667/ 

C5D6/. 83333333333/ 

C1D6/. 16666666667/ 

Cl  015/ . 66666666667E- 1 / 

C2D96/120. 42729108/ 

REPS/01  301000000000000/ 

RTOL/0 1 63 1 000000000000/ 

RK (  1  )/  .  1 6656066667 E  +  QO/ 

RK(  2 )/ . 53333333233E-01 / 

RK (  3 ) / . 2 1 333 33 3333 E +00/ 

RK (  4  )  / . 83233333233 E  +  OG/ 

RK (  5 )/. 26666666667 E+0 1 / 

RK (  6)/ . 25000000000E+01 / 

RK (  7 ) / .25781 2 5 0000 E+01 / 

RK  {  8  )  /  .  9  1  666*36666 7 E+01  / 

RK (  9 ) / . 66405250000E+01 / 

RK ( 10)/ .88541 666067 E+00/ 

RK( 1 1 )/ . 24000000000E+01 / 

RK (12)/. 80000000000 E+01 / 

RK( 13 )/ . 656045751 63 E+01/ 

RK ( 14 )/ . 3C555555556E+00/ 

RKf 1 5 ) / . 3 4509 30 392 2 E+00/ 

RK ( 16)/. 55086666567 E+00/ 

RKf 17)/ . 16533333333 E+01/ 

RKf 18)/ .94553823529E+00/ 

RKf 19J/.32400000000E+00/ 

RKf 20)/ .23273323529 E+00/ 

RK ( 21 )/ .2C354S51163E+01/ 

RKf 22)/ .69767441 860 E+01/ 

RKf 23)/ .56481 7981 46E+01/ 


0  /EK 2480 
D/EK2490 
D/EK2500 
D  J  E  K  2  5 ■  1  0 
D/EK2520 
D/EK2530 
D /EK  2540 
D;EK2550 
D/EK2560 
D7EK2570 
D7  EK  2580 
DVEK2590 
D/EK26C0 
D/EK2610 
DVEK2620 
D/EK2630 
0/EK2540 
D VEK  2550 
D;EK2660 
•DVEK2670 
DVEK2680 
D7  EK  2690 
Cy EK2700 
D  y  E  K  2  7  1  0 
D/EK272C 
D/EK2730 
DVEK2740 
D  y EK  2750 
D7EK2760 
DvEK2770 
D7EK2780 
Ov EK2790 
DVEK2800 
DVEK23  t  0 
0  7EK2820 
Cv/£K2830 
C7EK2840 
DVEK2850 
C  v EK2860 
DVEK2370 
CVEK2380 
CVEK2890 
D VEK  2900 
C  7  EK  29 1 0 
CVEK2920 
CVEK2930 
C7EK2940 
CVEK2950 
CVEK2960 
C  V  EK  2970 
CVEK2980 
CVEK2990 
C  VEK  3000 
DVEK3010 
C VEK  3020 
CVEK3030 
CVEK3C40 
CVEK3050 
CVEK3060 
DVEK3070 
C VEK3030 
CVEK3090 


oooooonoooooooooooooooooo  ooo  ooo  no  o  ooo 


DATA  RK ( 24 ) / . 1 3738 1 5676 1 E +00/ 

OATA  RK( 25)/ . 286302266 1 0E+00/ 

DATA  RKC26)/. 1 44 1 7355672 E+GO/ 

DATA  RK ( 27 ) / . 7 50000 0 CO 00 £ -01 / 

DATA  RK{ 28 ) / . 389928698756+00/ 

DATA  RK ( 29 ) / .31 944444444 E +00/ 

DATA  RK { 30 ) / , 135038363176+00/ 

DATA  RK ( 31  )/. 1 C783298827 E -  01  / 

DATA  RK<  32)/ . 698051948056-01/ 

DATA  RK ( 33 ) / . 62500000000 E -  02/ 

DATA  RK ( 34 ) / . SS6301 24777 E- 02/ 

DATA  RK { 35 ) / . 69444444444 E -  02/ 

DATA  RK { 36 ) / . 6 1 38 1 074 1 69E -  02/ 

DATA  RK ( 37 ) / . 68 1 8 1 8 1 8 1 32E -  01  / 

DATA  RK { 38 ) / . 1 0783298827E -  01  / 

DATA  RK ( 39 ) / . 69805 1 94805 E -  01  / 

BEGIN  INITIALIZATION,  PARAMETER 
CHECKING.  INTERRUPT  RE-ENTRIES 

IER  =  0 

ABORT  IF  IND  OUT  OF  RANGE  1  TO  6 
IF  (IND. LT.1 .OR. IND. GT. 6)  GO  TO  290 

CASES  -  INITIAL  ENTRY,  NORMAL 
RE-ENTRY.  INTERRUPT  RE-ENTRIES 

GO  TO  (5, 5, 40, 1^5. 265, 265)  .  IND 

CASE  1  -  INITIAL  ENTRY  (IND  . EO -  1 

OR  2)  ABORT  IF  N . GT . NW  OR  TOL.LE.O 

5  IF  (N.GT.NW.OR.TOL. LE.ZERO)  GO  TO  295 
IF  (IND. E0.2)  GO  TO  15 

INITIAL  ENTRY  WITHOUT  OPTIONS  (IND 
'.EQ.  1)  SET  C(1)  TO  C(9)  EQUAL  TO 
0 

DO  10  K= 1  .9 
C(K)  =  ZERO 
10  CONTINUE 
GO  TO  30 

SUMMARY  OF  THE  COMPONENTS  OF  THE 
COMMUNICATIONS  VECTOR 
PRESCRIBED  AT  THE  OPTION 
r-  OF  THE  USER 

CM)  ERROR  CONTROL  INDICATOR 
C(2)  FLOOR  VALUE 
C( 3  )  HMIN  SPECIFICATION 
C(4)  HSTART  SPECIFICATION 
C<5)  SCALE  SPECIFICATION 
C ( 5 )  HM AX  SPECIFICATION 
C( 7)  MAX  NO  OF  FCN  EVALS 
C(8)  INTERRUPT  NO  1 
C ( 9  )  INTERRUPT  NO  2. 

DETERMINED  BY  THE  PROGRAM 

CMO)  RRE3  (  RE  L  ROUNOOFF  ERROR  BND  ) 
C( 1 1 )  DWARF  (VERY  SMALL  MACH  NC ) 

C( 1 2  )  WEIGHTED  NORM  Y 
CM3)  HMIN 
C ( 1 4 )  HMAG 
C (  1 S  l  SCALE 
CM  6)  HMAX 
C (  1  7 )  XTRIAL 


D7EK3100 
DVEK31  10 
D/EK3120 
D./EK3130 
d;ek314o 
DVEK3150 
DVEK3160 
0;EK3170 
DVEK31SC 
D/EK3190 
D/EK3200 
D;EK3210 
DVEK3220 
DVEK3230 
D VEK3240 
DVEK3250 
DVEK3260 
DVEK3270 
DVEK3280 
DVEK3290 
DVEK3300 
D  V  E  K  3  3 1  0 
D  VEK  3320 
0 VEK 3330 
DVEK3340 
DVEK3350 
DVEK33SO 
DVEK337  0 
DVEK33S0 
D  V  EK  3390 
D  VEK  3400 
DVEK3410 
DYEK3420 
CVEK3430 
DYEK3440 
DVEK3450 
DVEK34S0 
DVEK3470 
C  v'  E  K  3  4  S  0 
D  VEK  3490 
CVEK3500 
C VEK 35 1 0 
uVEK  3520 
CVEK3530 
CVEK3540 
DVEK2550 
CVEK3560 
CVEK3570 
CVEK35SO 
CVEK3590 
CVEK3600 
C  VEK 361 0 
C  VEK  3620 
DVEK3630 
C VEK 3640 
CVEK3650 
DVEK3660 
DVEK3670 
CVEK3680 
D V  EK  3690 
CVEK37C0 
CVEK3710 


ooooooooo 


Appendix  D 


Material  Property  Data 

It  was  necessary  to  perforin  a  series  of  elastic 

modulus  determination  tests  to  characterize  this  adherent 

material.  Slight  variations  in  material  properties  can 

be  evident  in  molding  compounds  even  manufactured  by  the 

same  supplier.  In  a  separate,  extensive  study  concerning 

material  property  data,  Taggart  reported  the  elastic 

6 

modulus  of  SMC-25  to  be  2.1x10  psi  and  the  results  shown 
in  Table  4  are  in  close  agreement. 


Table  4 


Specimen  # 
SPEC1 
SPEC2 
SPEC3 


Modulus  (PSI) 
2 . 21xl06 
2.26xl06 
2.18xl06 


Plate  6  shows  a  typical  test  specimen  used  for  modulus 
determination.  The  data  for  these  tests  may  be  found 
on  the  following  pages. 
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Specimen 
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(!SM)  -o 


Microstrain 


M  icrostra  in 


Specimen 
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M  icrostrain 


Specimen 
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" Joggle-Lap"  Joint  Subject  to  Tension 


PLATE  4 
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PLATE  7 :  Photomicrographs  Showing  Relative  Fiber  Content 


